THE BULLETIN OF 
Mathematical 


BIOPHYSICS 


JUNE 1955 = 


a 


_ Analysis of Tracer Experiments in Nonconservative Steady-state Systems— 
HL E. Hart ie oe = = a ad - cd - a ad = - - - - - - e = 87 


__A Countercurrent Mechanism with Transverse Diffusion—George Karreman- 95 


_ Simultaneous Nonlinear Equations of Growth—W. J. Cunningham Rigsitee SE 


i Some Theorems in Topology and a Possible Biological Implication—N. Ra- 


Alka Mesias esis a cl Sa a Se eae TTT 
Nets on the Periodic Properties of Biochemical Systems Represented by 
Linear Equations—H. D. Landahl rr a 0/7 f 


B A Mathematical Model for the Temporal Pattern of a Population Structure, _ 
with Particular Reference to the Flour Beetle: Il. Competition between 
Ropecwereafandaht << = = ek en 2 181 


Diffusion and Simultaneous Chemical Reactions: |. A Method for Solving the 
Equations of Some Systems in Which a Fixed Concentration Exists at a 
Boundary—D. G. O’Sullivan- - - - - - - = = * = = = = = 141 


PRESS - CHICAGO 
- NUMBER 2 


: _ Some Order-Disorder Considerations in Living Systems—Harold J. Morowitz <>81 x sh 
' [ome ie ate 


The Bulletin of 
MATHEMATICAL BIOPHYSICS 


EpItTor: 


N. RASHEVSKY, UNIVERSITY OF CHICAGO 
Assoctate Eprrors: 
H. D. LANDAHL, UNIVERSITY OF CHICAGO 
ANATOL RAPOPORT, UNIVERSITY OF MICHIGAN 


The BULLETIN is devoted to publications of research in Mathematical 
Biology, as described on the inside back cover. 


THE BULLETIN is published by the University of Chicago at the University of 
Chicago Press, 5750 Ellis Avenue, Chicago 37, Illinois, quarterly, in March, June, 
September, December. {The subscription price is $10.00 per volume, the price of 
single copies is $3.00. Orders for service of less than a full volume will be charged at 
the single-copy rate. Postage is charged extra as follows: for Canada, 20 cents per 
volume (total, $10.20); for all other countries in the Postal Union, 50 cents per vol- 
ume (total, $10.50). {Patrons are requested to make all remittances payable to the 
University of Chicago Press in postal or express money orders or bank drafts. 


THE FOLLOWING is an authorized agent: 


For the British Commonwealth, except North America, India, and Australasia: 
The Cambridge University Press, Bentley House, 200 Euston Road, London, N.W. 1. 
Prices of yearly subscriptions and of single copies may be had on application. 


CLAIMS FOR MISSING NUMBERS should be made within the month following the 
regular month of publication, The publishers expect to supply missing numbers free 
~ only when losses have been sustained in transit, and when the reserve stock will 
permit. 


BUSINESS CORRESPONDENCE should be addressed to the University of Chicago 
Press, Chicago 37, Ill. 


COMMUNICATIONS FOR THE EDITOR and manuscripts should be addressed to N. 
Rashevsky, Editorial Office of The Bulletin of Mathematical Biophysics, 5741 Drexel 
Avenue, Chicago 37, Ill. a 

~NOTICE TO SUBSCRIBERS 
If you change your address, please notify us and your local postmaster immediately. 
[Copyright 1955 by the University of Chicago] 


Permission to reproduce for purely scientific or scholarly purposes any papel 


published in this journal will be given upon meget 


ee 


c 
a 


BULLETIN OF 
MATHEMATICAL BIOPHYSICS 
VOLUME 17, 1955 


SOME ORDER-DISORDER CONSIDERATIONS IN 
LIVING SYSTEMS 


HAROLD J. Morow1tz* 


LABORATORY OF TECHNICAL DEVELOPMENT, NATIONAL HEART INSTITUTE, NATIONAL IN- 
STITUTES OF HEALTH, PuBLIC HEALTH ‘SERVICE, U.S. DEPARTMENT OF HEALTH, 
EpucaTion AND WELFARE, BETHESDA, Mp. 


Order and disorder in biological systems are considered quantitatively in terms of infor- 
mation and entropy. After discussing the factors contributing to the information content of a 
living cell, a calculation is made of this parameter. The value for a typical bacterial cell is 
4.6 X 10!° bits. This value is compared with an experimental value of the heat of growth and 
entropy production of £. Coli. A discussion of methods of improving the calculation is also 
presented. 


While it has long been recognized that biological systems possess a 
high degree of order, it is difficult to define the degree of order in a 
quantitative fashion. This communication suggests a method of calculat- 
ing the degree of order and discusses certain consequences of this ap- 
proach. 

As a starting point let us consider a cell with a specified volume and 
atomic composition. It is physically possible for this system to exist in a 
very large number of configurations VV. Of these configurations a much 
smaller number, L, correspond to the cell being alive. If all states have 
equal a priori probability, then the quantity L/N gives the a priori prob- 
ability of the system constituting a living cell. If instead of assuming that 
all states have equal a priori probabilities, we consider p, as the a priori 
probability of the 7th state, then the probability of a living coaie uta yon, 


p, is given by 
OinPi 
abut (1): 
p aa =P; 


where 6;z, is unity when the 7th state is a living one and zero for all other 

states. The quantity p, L/N, or some function of this quantity represents 

a measure of the order of the system. 
The idea of representing order in terms of some function of the number 
* Present address: Sloane Physics Laboratory, Yale University, New Haven 11, Conn. 
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of states of the system suggests two familiar functions, entropy and infor- 
mation, as methods of making a quantitative evaluation. The entropy of 
formation AS and the information content of the cell, J, in terms of L 
and J are as follows: 
AS = k log L — k log N 
(2) 
l= loge NV — logs ds 


where AS is given by the familiar Boltzmann formulation (& being the 
Boltzmann constant) and J is the information content of the message 
that a given configuration corresponds to the cell being alive. 

The information concept has certain advantages in that no assumptions 
need to be made about equilibrium distributions of the states under 
consideration. For this reason information will be employed in subsequent 
calculations. 

Before proceeding to a more detailed treatment, we shall consider cer- 
tain experimental results which simplify the procedure. The first experi- 
mental finding is that certain unicellular organisms such as bacteria are 
able to withstand almost complete drying. Routine lyophilization pro- 
cedures remove up to 95% of cellular water without inactivating the 
organisms. Since adding water to the dried cells is a random process 
which does not add information to the system, we may conclude that the 
basic information unit of living cells is in the nonaqueous portion. The 
second experimental finding is the ability of certain organisms to survive 
temperatures near the absolute of zero. Recently the author and Mr. 
William Hanson were able to subject dry B. subtilis spores to tempera- 
tures of approximately 1.3° K for a period of 53 hours. A high percentage 
of the spores survived this treatment, indicating that the basic informa- 
tion unit is in the structure quite independent of the vibrational and 
rotational state of the system. Thus for computational purposes we may 
consider the nonaqueous portion of the cell at absolute zero as constitut- 
ing the basic information unit necessary to specify a living system. 

To evaluate J in equation (2) we shall consider that log: L is negligible 
in comparison to logs V. This is equivalent to assuming that only one state 
corresponds to the system being alive. Since L is very small compared to 
N this is likely to cause little error in J, although the procedure is nonethe- 
less somewhat unsatisfactory since the quantity L would be a very inter- 
esting one to be able to calculate. In any case if one makes the above as- 
sumption, equation (2) takes the simple form 


I=logN. (3) 
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The quantity NV, it will be remembered, is the number of structural 
states corresponding to all possible positions and bonds, the system being 
constrained by atomic composition and volume. 

Consider a system containing 1, of the ith type of atom with a valence 
of a;. The total number of atoms in the system M is 


Mi p> NG. (4) 
The total number of covalent bonds which may be formed is 
ee 
Ca ain; . (5) 


We shall proceed to estimate NV, the total number of possible configura- 
tions of the system. 

Since we are considering the system at absolute zero, we need consider 
only the coordinates and chemical bonds of each atom. If we further as- 
sume that each atom occupies approximately the same volume, we may 
consider the cell as being made up of a number of boxes (unit atomic 
volumes) with the first approximation of the number of states, V1 being 
given by the number of permutations of the M@ atoms in M boxes: 


A state so defined does not consider the number of bonding states 
possible within a given configurational state. The following illustrative 
diagram shows a number of bond states which may exist for a given con- 
figurational state with this same set of atoms: 


H H-O H-H 0O H H-O 
vib | 

O=C H aoe H ‘eae 

alee H-—C-H H O-H. 


If an atom of the ith type can distribute its bonds among its neighbors 
in B; independent ways, then the number of bond states for a given con- 
figuration is less than B, where B is defined by equation (7): 


B=] | Bi: (7) 


The number of states V is then smaller than 1; B and ) 
log N<log NitlogB. (8) 
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Using Stirling’s approximation we get 


log NM log M— do 1: log nit >) ms log B;. (9) 


If we consider each atom as having six nearest neighbors, we get By = 
6, Bo = 21, Bc = 120, By = 56, Bp = 252, and Bs = 21, where the sub- 
script denotes the chemical symbol of the species of atom under con- 
sideration. 

To obtain a numerical value of the order of magnitude of J, we perform 
the calculation on a structure the size of a small bacterium, the mass 
of the nonaqueous portion being 10-18 grams. The atomic composition of 
a typical bacterial cell of 10-"% grams dry weight is given in Table I. 


TABLE I* 
NN 6ADAG O 4.7 X 108 
C24 16 Peo 10 
H PAZ SAY" Se sre 10? 


* Atomic composition of a single bacterium (10718 grams 
dry weight) in terms of numbers of atoms. 


Using these values we arrive at a value of J of 4.6 X 10% bits. While the 
assumptions seem somewhat arbitrary it should be realized that the 
assumed values may cover a wide range without causing an order of mag- 
nitude change in J. 

The assumption of equal a priori probability of all states poses the most 
serious difficulty in the preceding analysis. 

Following equation (1), the information content of the system under 
consideration is given by equation (10), when no assumption is made 
regarding equal a priori probability 


[== D3 (pi loge pi — bixP: logs pi) « (10) 


For an equilibrium distribution of states p; is given by the following 
expression * 


Se -B/ kT . (11) 


All energies may be taken relative to a zero value for the system with 
no bonding. The #;, are therefore negative quantities. The energies under 
consideration are equivalent to free energies of formation at absolute zero. 
The difficulty in making calculations based on equations (10) and (11) is 
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our lack of knowledge regarding the energy distributions of the various 
states, both living and non-living. 

Two pathways seem available for achieving a better estimate of the 
effect of the probability distribution of different states; one involves a 
detailed theoretical analysis of energy distribution in different states of the 
system while the second deals with experimental attempts to get the free 
energy of formation of such structures as dry spores or other materials 
that can be handled by calorimetric techniques. 

Returning to equation (2) we note the relation between information 
and the Boltzmann definition of entropy 


AS = —kI log, 2. (12) 


Consider now the growth of bacteria in an isothermal calorimeter. For 
each bacterium produced there is a local increase of information J. In 
accordance with the second law of thermodynamics there must be an 
increase of entropy of at least AS somewhere else in the system. If heat 
is transferred slowly to a sink at temperature 7, then the entropy change 
will be in the order of AQ/T. Thus an estimate of the maximum informa- 
tion content may be obtained from 


a= — BF log. 2. (13) 


Values of AQ have been measured for bacteria by S. Bayne-Jones and 
R. S. Rhees (1929). The AQ associated with the production of one organ- 
ism varies during the growth cycle, but we shall here consider the mini- 
mum value obtained. For B. coli a value of about 4 X 10-1? gm. calories/ 
organism is gotten at 37.5°C from calculations based on the data of 
Bayne-Jones and Rhees. The minimum value of AQ based on the previous- 
ly obtained value of the information content (corrected to the mass of a 
B. coli, which is approximately 4 X 10- gms. dry weight) is about 
1.2 X 10-” gm. calories/organism. The heat actually given off is about 
three times that required to account for the decrease of entropy. Or, 
stating the result alternatively, the information content cannot be more 
than three times that calculated by the preceding theory. 

_ However, since the bacterium is not 100 per cent efficient in ordering 

‘the constituent components, we may interpret the result to read that less 
than 30 per cent of the heat given off by the growing bacteria is the result 
of local decrease of entropy, while 70 per cent is the result of other 
processes such as motility, maintenance of concentration gradients, and 


various metabolic processes. 
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Various factors in calculation of the information content tend to give 
a maximum possible value of this quantity. Among these factors are 
assuming equal a priori probability, taking the maximum value of B, 
and neglecting log L in comparison with log V. The order of magnitude 
of the calculated value of AQ in comparison with experimental values, 
however, indicates that the computed value of J is in a reasonable range. 
Further work is necessary to arrive at a more detailed value of the infor- 
mation content. 


LITERATURE 
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A method is developed for finding the transfer and localization rates and the volumes of 
N compartment steady-state biological systems from experimental results. It is shown that a 
complete solution for certain systems in which the rates and volumes remain constant and in 
which there is access to all compartments can be obtained by using a single radioactive tracer. 
The information obtainable from experiments wherein some compartments are not accessible is 
analyzed for mammillary and catenary systems. Conservative systems are handled as special 
cases in which the localization is zero while anisotropic membranes separating compartments 
are shown to introduce no additional mathematical difficulty whenever all compartments are 
accessible. The limitations on the use of this method of multi-compartment tracer analysis 
are briefly discussed. 


Introduction. Many studies of compartment volumes and exchange 
rates of physiological and non-physiological substances in biological 
systems have appeared in recent years (Gellhorn, Merrell, and Rankin, 
1944; Flexner, Cowie, and Vosburgh, 1948; Solomon, 1949; Sheppard and 
Householder, 1951). 

In general, these studies have been concerned with steady-state or 
quasi-steady-state systems in which the labeled material is either a non- 
localized or very slowly localized dye or radioactively tagged substance. 
It will be the primary purpose of the present work to analyze tracer 
experiments with localizable substances in NV compartment systems of 
which M compartments are available for measurements, where all volumes 
and rates remain constant and under conditions in which experimental 
errors are small. 

In Section I, the complete solution of an N compartment non- 
conservative steady-state system is obtained assuming access to all 
compartments. In Section IT, the special cases of mammillary, catenary, 
conservative, and anisotropic systems are treated. In conclusion, the limi- 
tations of this and other tracer methods are briefly discussed. 


* Research supported by the Atomic Energy Commission, Contract AT (30-1)-1551. 
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TABLE OF NOTATION 

x“; = concentration of freely diffusible labeled material in the ith com- 
partment.* 

x;1 = concentration of freely diffusible labeled material in the 7th com- 
partment at time ¢;. 

V; = volume of ith compartment. 

r,; = rate of fluid transfer from the ith to the jth compartment. 

—r,; = the rate of localization or precipitation of the tagged material 
within the 7th compartment. 

—a;; = ri;/V; = rate of fluid transfer from the ith to the jth compart- 
ment per unit volume in the 7th compartment. 


—A;4; = «(> rs)/ V; = rate of transfer per unit volume in the 7th 


compartment of labeled material from the 7th to all other compartments 
plus rate of localization per unit volume within the ith compartment. 

Axm = coefficient of et in the expansion of xx in terms of exponential 
functions. 


t 
Vi = uit es ij fe x,(t)dt = concentration of labeled material present 
7 0 
in all forms in compartment 7 at time #,. 


I. Measurements in all N compartments 
Let x; be the concentration of freely diffusible labeled material in the 
ith compartment, V; the volume of the ith compartment, 7;; the rate of 
fluid transfer from the ith to the jth compartment, —r;,«; the total rate 
of localization or precipitation of the tagged material in the ith compart- 


ment; then, assuming 7;; = 7;:, it follows from the equation of continuity 
that 


U 


ie 
t= >> 
7=1 


Since the labeled material is present only in tracer quantities, 7;; and 
rj are independent of x; and %;. It will be assumed that 7;;, 7;;, and V; 


N ? 
rig (Xj — 4%) + rex, where bs. indicatesi#~ 7. (1) 
7 


* When studying a physiological material, unlabeled as well as tagged molecules will always 
be present and it may then be convenient to consider x; as the specific activity in w¢/mg. of 
those molecules of the material under investigation which are present in freely diffusible form 
in the ith compartment; V; would then be the quantity in mg. of such non-localized molecules 
in the ith compartment and 7;; would be the rate in mg./unit time at which the material is 
transported from the 7th to the jth compartment. 


} For high molecular weight substances which do not readily diffuse through compartment 
membranes, the 7;; are in the nature of clearance rates rather than fluid transfer rates. 
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are constant. Let 


N 
fea—V; >) aij, (2) 
7=1 
hei) a V 5a; , AG ) (3) 
and equation (1) becomes 
ts+ > ai3X; = 0. (4) 
7 


If simultaneous measurements of the concentrations in all compart- 


ments can be made repeatedly, equation (4) can be expressed in terms of 
experimentally determined quantities: 


Lip > a; j%j,= 0, (5) 
Fi 
where 4; is the measured concentration of tracer in the ith compartment 


at time ¢;, and #;; may be determined graphically or algebraically. 
Solving equation (5) for the transfer rates per unit volume, 


Nye - » H{-1,1 — Li Nj4ij1+ + + XN1 
X12 Me woy cae Aen Xe 

aij= F Spi) Ge (6) 
Xin ny = a ee Gee 


Since the tracer material must be conserved, 
N ty 
ae 1 Viti—7 f (t) 


and K is the amount of tracer originally introduced. From equations (2) 
and (7) it follows that: 


Yu--- Yi-1,1 1 Viti1l-+-YN1 


712 ViH—1,2 1 Yi41,2 YNn2 
Wes TOE + \¥ (8) 


vin Yi-1n 1 Yit1Nn aa | 
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where 
t] t] 
y= eat Do ai; f x; (t) dt and ip x; (t) dt (9) 


are evaluated numerically. 

Equations (6), (8), (3), and (2) can therefore be used to determine the 
transfer rates per unit volume a;;, the volumes V;, the transfer rates 7;;, 
and the localization rates 7;; for an N compartment system. All that is 
required is a single radioactive isotope or other label, and access to all 
compartments so that a sufficient number of-simultaneous measurements 
of the concentrations can be performed. 

A mathematically equivalent but more indirect method can be used in 
which the measurements need not be simultaneous (although repeated 
access to all compartments is still required). From the theory of simul- 
taneous first order linear differential equations (Ince, 1926) it follows that: 


d 
Qi ea a2 --.- Qn 
d 
G21 a9 + ira don 
“,=0, Pitzer pee (L@) 
a a + g | 
N1 N2 ann ai | 
Rewriting (10) 
) r > 
Be Cee? Ca 5 th cn a Creat + Cyr, = 0, (11) 


where «{” is the Pth time derivative of x; and C, are known functions of 


a;;. Furthermore 
we= >> Anve™** (12) 
™ 
where dy are the roots of the secular equation 
AV +CiN 14 CN 2+. Cy a Cy RO. (13) 


Substituting expansions (12) in (4) and setting the coefficients of 
e\“' = () one obtains the NV? relations 


lu Am+ > aj;4ie=0, i,j, M=1...N, (14) 
7 ? 


* It will be assumed in this treatment that there are no multiple roots of 2. However, not 
all exponential terms need appear in each of the compartments (see Section IIb). 
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therefore, 
eet rt, — Ni an Atte el 


Aj Age Aj-12 = 2A 72 Aji Ano 
ajj= + |A (1:5)) 


A 1N A 2N A j-1,N = AwA tN A j+1,N A NN 


The volumes V; can still be determined by using (8) with 7; redefined: 


= Aye] ( eu = 1) : 
Yit = — Agee = Ss ai; Ady leah —t : (16) 
; 


M Me 


The numerical values of A xy, and Ay in (12) can be obtained graphically 
by well known methods whenever the roots of (13) are moderately spaced 
and the measurements are carried out over a sufficiently long period of 
time so that eventually all but one exponential term are observed to 
become negligible.* Since the values are determined graphically, a large 
number of experimental points can be used to help ascertain the slopes 
and intercepts of the various straight lines. The resulting values of A xyz 
and )y, are in the nature of average values and the effect of small random 
experimental errors is somewhat reduced by this method. The over-all 
consistency of data can be checked immediately since the graphically 
determined values of \, should be the same for all compartments. Finally, 
of course, measurements in different compartments need not be made 
simultaneously. 

Both methods outlined above are essentially equivalent and the ex- 
perimental conditions necessarily determine the procedure of choice. 

It can be shown that lack of access to a compartment limits mathe- 
matical analysis, but if enough supplementary facts are known, such as 
some of the volumes or transfer rates, then a complete or partial descrip- 
tion of the system may still be possible. This is particularly true for con- 
strained systems in which many of the transfer rates vanish and such 
special systems will now be used as illustrations. 


II. Special cases 


a) Catenary systems 
The condition that 7,; = 0 for 7 #i — 1, 7,¢-+ 1 defines a catenary 


* When two or more roots da, 4 are sufficiently closely spaced, accurate determination of 
a; and V; may be impossible because of the indeterminacy in Axa and Axg. 
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system. Equation (5) becomes: 


deat Oita + eintir + Giyavtii = 0, %1=1,2,3. (17) 


The transfer rates per unit volume in the 7th compartment (ie., ais, 
@;,:41, and a,,;1) can therefore be found from equation (17) using repeated 
simultaneous measurements in the 7 — 1, 7, andz + 1 compartments. 

b) Mammillary systems 

A mammillary system in which compartment 1 is the central pool is 
defined by the requirement that 7;; be non-vanishing only if 7 = 1 or 
4 = 1 ori = j. Equations (14) become: 


(Au + @xx) Axr+@niAm=0, KA1. (18) 


If compartments 1 and K are accessible, Axy, Aim, and Ay will be 
measurable. Therefore, axx and @x can be determined, a result previously 
obtained (Sheppard and Householder, Joc. cit.). 

That the number of exponential terms in the expansion (12) for any 
one compartment may be less than the number of compartments actually 
present is readily demonstrated for a mammillary system. Solving for 
Axv in equations (18), substituting the results in (15), and simplifying the 
resulting determinants by the repeated application of standard methods 
one obtains: 

Oitéx) eter)... Awtax) 


axk= 
Oxi (d2— Ax) (d3—ax)...(@x-1— @x) (4x41—Gx)...(@y—Ax) ’ 


(19) 


where a; is written instead of a;;. 

If @1t = xx, t ¥ 1, k, the denominator of equation (19) vanishes and 
since dix is not infinite, it follows that \y, = —daxx for at least one value 
(M = P). From (18) it follows that since Axp is also not infinite, Aip 
must vanish. Since Aip = 0, (18) further implies that (Ap + a;,)Aip = 0, 
whence A;p = 0 for all values of 7 such that a; ¥ dxx. If all Aip (i ¥ &, 6) 
vanish, then measurements in the central pool and all other compart- 
ments 7 for which a;; # @xx would give rise to at most NV — 1 exponential 
terms in each of the expansions (12) and would therefore superficially 
suggest that the system had only N — 1 compartments. Since 


tte treet? Ri 
= Se ca ee Fo 
it follows that a;:= axx does not imply that 
Pies fe Cit oo hey 
Vs K Hs V; K. (21) 
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Both compartments (i and k) can be quite different and yet one less 
exponential term will appear in the central compartment. 

c) Conservative systems 

Such systems may be defined by the condition that all 7,; = 0. If the 
volumes are known initially then access to only V — 1 compartments is 
required inasmuch as equations (7) can be used to complete the solution. 

d) Anisotropic membranes 

Equation (1) assumes that there are free flow and diffusion between 
compartments. If the conditions are such that there is a preferred direc- 
tion in fluid flow between compartments then the volumes are likely to 
vary (especially in non-biological systems) and steady-state methods 
may no longer apply. However, if the volumes remain constant, steady- 
state methods can be applied and it follows that: 


/ 
Vit;= > (75X53 — rigs) + rita. (22) 
- 


Setting r{, = rus — we '(r3s — raj), equation (22) can be rewritten as 


7 
/ 
ai >> T 54 (4%; — X;) + risks - (23) 
7 


Equation (23) is formally equivalent to equation (1) and the methods 
previously developed can be used to find 7;:, 7:;, and V; for all 2 and 7 
whenever there is access to all compartments. * 


ITI. Discussion 
In the above treatment, only measurements of freely diffusible tracer 
material are employed. It has been pointed out that measurements of 
the localized material may on occasion also be used to advantage.} Thus: 
M (t) 


: 24 
fxd ( ) 
0 


where M(t) is the amount localized in the 7th compartment. 

From equations (6), (15), and (8) it is readily seen that if the value of 
any determinant is small compared to individual elements of the de- 
terminant, small percentage errors in the elements may result in a large 


gS 


* The equation of conservation of tracer material becomes 


a { Viti rh f 'xdt} =K, 


u 


+ Professor H. D. Landahl. Private communication. 
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error for the calculated rates and volumes. This uncertainty in the calcu- 
lated rates and volumes will also tend to increase rapidly as NV increases. 
Furthermore, since no biological system is absolutely static the various 
rates and volumes may themselves be changing. For such quasi-steady- 
state systems or systems in which repeated measurements are not feasible 
the use of multiple tracers is indicated (Sheppard and Householder, loc. 
cit.). 

Finally, it should be noted that any purely tracer method gives little 
information on the non-linear properties of a system. If a tracer X is 
added to a steady-state system, the resulting differential equations in X 
will usually, if not always, be linear with constant coefficients and conse- 
quently give little information on the non-linear equations which deter- 
mine the steady-state concentrations of the stable or unlabeled material. 
The use of modified tracer methods in studying non-linear systems would 
seem worthy of investigation. 


I am grateful to Dr. Daniel Laszlo for his stimulating encouragement 
and support of this work. It is a pleasure to thank Dr. Robert Hatcher 
and Mr. Wallace Gold for their careful review of the theory and calcula- 
tions. Thanks are due to Mr. Mones Berman, Mr. Edward Siegel, Dr. 
Paul Weisz, and Dr. Aaron Yalow for helpful discussions. 
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A simple treatment is given for the stationary state distribution of a substance dissolved 
in two liquids flowing in opposite directions in two tubes with a common wall that is permeable 
to the substance. The transverse diffusion is taken into account. It is shown how in the case of 
sufficiently fast diffusion the permeability constant of the wall can be obtained from the con- 
centrations of the substance at the ends of the tubes, the velocities in the tubes, and the 
length of the tubes. 


An interesting countercurrent mechanism occurs in the gills of fishes 
forming a highly efficient diffusion apparatus for the exchange of oxygen 
between the blood in the capillaries of the lamellae of the gill-filaments of 
the gill-arches and the water between those lamellae (van Dam, 1938, - 
pp. 78-79). It is the purpose of the present work to give a simple theo- 
retical treatment of such a countercurrent mechanism, taking into account 
the transverse diffusion of the substance. 

To simplify matters it will be assumed that both the capillary and the 
adjacent water channel are straight tubes with rectangular cross-sections 
F, and F2 (Fig. 1) and that the capillary is in contract through a perme- 
able membrane with the water channel only at one side PORS and similar- 
ly for the latter. Furthermore we will assume that the permeability of the 
wall PORS to the substance is the same everywhere along the tubes. We 
will use the approximation method of N. Rashevsky (1948). The velocity 
of the water is denoted by 2, and that of the blood in the capillary by ». We 
then have in Figure 1 7, > 0 and » < 0. The height of the water channel 
is denoted by ay, that of the capillary by a2, whereas 6 represents the com- 
mon width of both tubes. We will subdivide both tubes, with cross-sec- 
tions ab and ab respectively (Fig. 1), into layer of thickness Ax, The 
average concentration of the substance in the water in such a layer is 
represented by %(«, é) and similarly that in the blood by (x, ¢), a being 
the coordinate in the direction of the axis of the tubes and ¢ the time. 

Considering only convection in the direction of x and a diffusion cross- 
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wise, we see that a layer of length Ax of the water channel loses an amount 
of substance due to convection given by 


0G, 
a,b 01{ G1 (a, t) — t (w-+Aa, t) } = —abn 2 Ax, 


whereas due to the flow in the direction perpendicular to the wall PQRS 
it loses 


bAx D, 21 = 2 bAxD, 


According to the principle of the material balance the sum of these two 
amounts each with a negative sign must be equal to the average rate of 


L 


FicureE 1 
change in time of the amount of substance @,a;bAx in the layer considered. 
This consideration leads us after some simplification to our first equation: 


(a1— cf) =a $4, (1) 


Ox ay, 


Similarly we obtain for the average concentration % of the substance in 
the capillary: . 


oop) 2D, 


0¢ 
— av barges 1 cr sor ds (22— G2) = a, —— ae) 


art - (2) 


In these equations D; and D denote the diffusion coefficients of the aun 
stance in the water and in the blood respectively, whereas 2; and é 
represent the average concentrations at the wall PORS of the substance 
in the water and in the blood respectively. The continuity of the flux in 
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the transverse direction is expressed by: 


, = 
Di = hi th) = DA (3) 


in which represents the permeability of the wall PQRS to the substance. 
In the stationary state we have: 


Jun Oc 
baleeeeas ti (4) 
Introducing (4) into (1) and (2) we obtain: 
LCi ee. 
@, 0} oe ted) 4 (5) 
and 
dé PAY) 
— 2% Fo + (eh — ta) = 0 (6) 


since in the stationary state and @ are functions of x only. From (3), 
(5), and (6) we find: 
—_ = gq. ——_ ee 


ae Oe 
Se (8) 
Integration of (7) yields: 
Ci = 0 Ce + A (9) 


in which A is an integration constant determined by the boundary condi- 
tion that ¢ = Gio and & = Gp for x = 0. We then find: 


A = Cy — alan. (10) 

Introducing (10) into (9) and rearranging we obtain: 
61 — ACo = C10 — 2060 (Geb) 

or, together with (8), the physically obvious relation 
01 0161+ Gy V2 Eq = 1 V1 C19 + Aq V2 E29 « (12) 


This relation between the concentrations @, and @ is a SEN condition 
for such a mechanism considered here. 
Us . 
Solution of the two equations (3) for ¢ and @, yields: 


_, _ Di (a2h + 2 Ds) ér+ ah Does (13) 
1 (Di a2-+ Dads) h-+ 2D;Dz 
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and 
Dz (a;h +2 Dy) 2+ ah Dex 


(Dia.+ D2a;) h+2D,Dz2 (14) 


=f 
62 = 


Solving (9) for @ and introducing the result into (13) and then the ex- 
pression so obtained for @; into (5) we have with (8) and (10): 


di 2D, Doh «nieve 
da: a0i| (Diack Dig) hy 2D Del ae 
(15) 
a 2 Di Deh (E19 — a E20) 
a, { (Dya2+ Day) h+ 2D; D2} a° 
Solution of (15) gives: 
= Ce (16) 
in which 
h —1 2D,D. 
Oa a ag ee (17) 
ay a (Dia2+ D2a;) h+2D,Dz 
or, using (8), 
= h aut Gg V2 2D,Dz (18) 


Q, 0; G2 V2 (Dia2+ Dea;) h+2D, Ds; 


and C is an integration constant determined by & = Go for x = 0. From 
this condition we find: 


C= (ei — é20) ; (19) 
or, with (8), 

_ G2 V2 Ps - 

“erica E20) « (20) 


Introduction of (20), (18), and (8) into (16) and the result into (11) and 
solution of the equation so obtained for @ gives the solution of our problem: 


z __ 4, 0, E19 + Ae V2 Ea9 G2 Vg (a Zon) 
1. = tt _ (ir -— 
211+ G2 2 2101+ G2 22 Ds 3 
(21) 
Ne e—2h(a,0, +4,0,) D, D,z/a,0,0,0,(D,agh+D,a,h+2D,D,) 
and 
ee 1 01 C10 + G2 V2 E20 Qi 21 (z Zoo) 
=> — ———__ (i — Ee 
2101 -F 2 Ve 01+ G22,” z= (22) 


>" 4 e—2h(ayv 1+4,v,) D,D,z/a,v 19%; (D,a,h+D,a,h+2D 1D) 
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If we assume to > tp» and a2, > —dy%», which is the case in the gilis of 
fishes, we see from (21) and (22) that both % and @ decrease as x increases 
from « = Q. Because the exponent of the exponentials in (21) and (22) 
is positive for x > 0 we see that if « is large enough and @ become nega- 
tive. Then since for that value of x for which @ = 0 we have (dé)/(dx) < 
0, therefore, we see that if the tubes are sufficiently long for this to happen 
we have to solve equations (5) and (6) under different boundary condi- 
tions, namely, if @ = 0 for x = x* then also (dz)/(dx) = 0 for x = x* 
(Rashevsky, Joc. cit., p. 29). However that would only yield the uninter- 
esting solution % = 0 and @ = 0 for all x. If the concentrations at the 
ends of both tubes are greater than zero as is, e.g., the case in the gills, 
formulae (21) and (22) hold for the whole length of the tubes. We then 
obtain by subtraction of (22) from (21): 


@:-—02= (E10 — C20) e—2%h(a,0,+4,0,) D,D,2/a,0,4,0,(D,a,h+D,a,h+2D,D,) (23) 


so that the difference in concentration of the substance in both tubes is 
an exponentially increasing function of x if to > to, v2 < 0 but if also 
202 — a0, > 0. 

If in the latter case the length of the tubes is / we find by introducing 
@ = &, and x = / in (21) and some rearrangement: 


1 l 
a2 ay 1 P a2 oe a, Vj (24) 
2Dit2Di' hk 2 V2 ( C19 — Ce0) ; 


log 1 01 (C11 — C10) + G2 V2 ( E11 — C29) 

To obtain the permeability # and the diffusion coefficients D; and D; 
we need two more relations between h, D:, and D:. The two equations 
(3) are two such relations, if @; and @, are known together with cor- 
responding @ and %. However to measure the former concentrations will 
of course only be possible if the tubes are sufficiently wide, therefore, it 
will be impractical in most physiological cases. If in such cases the dif- 
fusion is sufficiently fast, so that 


atees A 
apes 68 aD, Si 
(24) reduces to: 
= Lol Sle ot Re (25) 


2 V2 (E10 — E20) 


log 4 0; (E31 — E40) + 4202 (C11 — C20) 
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From (25) the permeability of the wall to the substance can be obtained 
from the dimensions of the tubes, the velocities in the tubes, and the 
concentrations of the substance at the ends of the tubes. 


The author wishes to thank Dr. L. van Dam for a discussion which has 
led to this work. 
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The simultaneous equations 


dx : 
a ON 
dy _ a, 
Gin Ele W@19 


have been suggested to describe competition between two species. Information about solutions 
for these equations is given as stability diagrams related to coefficients in the equations. Ex- 
amples of solutions are shown as found with an analog computer. 


Several pairs of simultaneous nonlinear differential equations have been 
proposed as describing the competition between two species of organisms. 
One of the simplest sets of equations studied in detail (Volterra, 1931) is 


dx a. 
—_ = 7 — _ 1 
dy Oy 

t Ry 


Here x and y are the populations of the two species and ¢ is time. Coeffi- 
cients az, dy, Bz, By, kz and ky are real constants. These coefficients deter- 
mine the nature of solutions for the equations. Among the possibilities are 
that one species disappears leaving the other, or that periodic oscillations 


occur in the populations. ' 
Another pair of equations (Gause and Witt, 1935) is 


yay 9 — Fit) 9. (4) 


The extra terms in these equations introduce a sort of damping effect at 


the expense of a more complicated solution. 
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A still more general pair of equations studied only briefly in the past 
(Hutchinson, 1947) is 


t=F [ke—¥— fe(9) 12, (5) 
y= a lky—9— Sula) 19. (6) 


Functions f,(y) and f,(x) are continuous, single-valued, and differen- 
tiable, and both functions vanish for zero argument. These functions 
might be polynomials such as 


fe(y) = Boy + yey" (7) 
tu) = Byx + yt" ) (8) 


where @z, By, Yz, and y, are real constants. Additional terms of higher de- 
gree might be present also. These functions provide the coupling between 
variables x and y. 

While solutions for x and y as explicit functions of ¢ appear to be quite 
difficult to find, considerable insight into the properties of these solutions 
can be found by studying the relation between x and y as ¢ varies im- 
plicitly. Such a study can be carried out by considering the equation ob- 
tained as the ratio of equations (5) and (6), 


of Ry — ae) 
dy _ By! y se Io) 
" Blke—@— fa(y)] 


Singular points (Minorsky, 1947), or points of equilibrium, exist where 
both numerator and denominator vanish simultaneously. Values of « and 
y at a singularity are designated as x, and y,. There are four groups of 
singularities, as follows, 


1.%=0,y=0, 
2. Xe = ke, Ye = 0, 
3.0, = 00S) = hy 
4, t= Re — falye)s Ye = Ry — fyl%s) « 


Each of the first three groups contains one singularity, the number of 
singularities in the fourth group depends upon the properties of functions 
f, and fy. 
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The nature of solutions near a singularity can be explored by replacing 
« with (wv. + u) and y with (y, + v), where wu and » are small quantities. 
At the same time, the coupling functions can be expanded as 


fe(Ys + 2) = folve) + felys)v +... 


fy (ts ate u) = fy(&a) + fi(xs)u + ae 
with 


4 acid 
eye) ~ Gy V2 (9) | ’ 
f(a) =£ fy (a1, 


If these substitutions are made, and only linear terms in u and 7 are re- 
tained, equation (9) takes the form 


dv _ Au+Buy 
du Cut Dv’ (10) 
where 
een CAN (11) 
B= 7" [hy — 2ye— fy (%) 1, (12) 
C =F [he Qt fe (ye), (13) 
D= —= w.ft (9s) . (14) 


The nature of solutions for equation (10), and thus of solutions for equa- 
tion (9) near a singularity, depends upon the characteristic roots 


(i, Ae) = 2{(B + C)  [(B+ C)? + 4(AD — BC)}}. (15) 


Real roots of opposite algebraic sign give a saddle point, which is always 
unstable in the sense that solutions ultimately grow without bound. Real 
roots of like sign give a nodal point, unstable only for positive roots. Com- 
plex conjugate roots give a focal point, unstable only if the real part of 
the roots is positive. A focal point represents an oscillatory solution, 
saddle and nodal points are nonoscillatory. 

Evaluation of the roots near each singularity gives the following results. 

1. Case x, = 0, ys = 0. Only coefficients a, and a, are important, and 


the stability is summarized in Figure 1. 
2, Case x, = kz, ve = 0. Evaluation gives B = (a,/ky)Y where Y = 
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ky — f,(Re), so that this quantity Y is important in determining stability. 
While Y can readily be found from its defining relation, a graphical inter- 
pretation may be helpful. The curve y = k, — f,(«), plotted in Figure 2, 
is the locus of values of « and y making the bracket in the numerator of 
equation (9) vanish. Any solution curve for equation (9) crosses this locus 
with horizontal slope; thus the locus is the isocline for zero slope, dy/dx = 0 
or y = 0. The ordinate of the isocline at x = k, is the quantity Y. The 
resulting stability is summarized in Figure 1, where now coordinate Q 


Va en ted 
N SS Ny N 7 
‘ ee N SS \ N Q Ya Ya 7 2 
XN \ \ \ \ 7 7 4 ce La 
S ~ N \ ¢ 7 we a 4 
\ S x Si 7 
Ss wesaddilesras) unstable 
Sah Ss 47s mode “ 
Ny \ ~ 7 vA 


\ 
a a aS ~ X N 


Go 7 3 
/? *estables SS 


N SS ‘ 
S Sataiellite SSS 
~ = \ 


\ Sy 2S ES 


SC eniGden ese ees 


Ficure 1. Stability of solutions near a singularity of any of the first three groups. Co- 
ordinates are: Group 1, P = az, Q = ay; Group 2, P =—az, Q = ayY/ky; Group 3, P =—ay, 
Q = azX/kz. 


x 
df 4 =O 


Ficure 2. Graphical determination of quantities X, Y, WM , and N, and location of singu- 
arities by means of the isoclines for y = 0 and ¢ = 0, 
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involves quantity Y. Among other things this diagram indicates that the 
algebraic sign of Y helps determine the sign of Q and thus the nature of 
solutions. 

3. Case x; = 0, y, = ky. This case is analogous to that just preceding. 
Evaluation gives C = (a,/kz)X where X = k, — f.(k,). The quantity X 
has the graphical interpretation also shown in Figure 2. The curve 
x =k, — f,(y) is the isocline of infinite slope, dy/dx = » ort =0. A 
solution curve for equation (9) crosses this line with vertical slope. The 
abscissa at y = ky is quantity X. The stability is summarized in Figure 1, 
where Q now involves quantity X. 

4. Case %. = kz — falc), Vs = ky — fy(xs). The number of singularities 
determined by these relations depends upon functions f, and fy. With non- 
linear functions, determination of «, and y, may be tedious, requiring solu- 
tion for the roots of an algebraic equation of high degree. Probably a 
simpler and more informative approach is through a graphical construc- 
tion. The singularities can be located at the intersections of the isoclines 
for ¢ = 0 and y = 0, as in Figure 2. Quadratic functions, as in equations 
(7) and (8), yield parabolic isoclines. The number of singularities produced 
by intersections of these two parabolas can be any integer from zero to 
four. Two such singularities should occur in Figure 2, although the figure 
is too small to show the one in the third quadrant. There are four inter- 
sections in the examples of Figures 4 and 5. More complicated coupling 
functions could yield more than four singularities. 

The values of quantities A, B, C, and D for this fourth group of singu- 
larities are most easily expressed in terms of additional abbreviations, 
R = a,%,/ke, S = ayy./ky, G = R/S, M = fy(x.), and N = f(y). Calcu- 
lation of these quantities requires that the singular point be known. 
Values of M and N can be found, if desired, from the slopes of the isoclines 
at the singularity, as in Figure 2. A consideration of equation (15) with 
these quantities substituted allows the following conclusions regarding the 
types of solutions. 


a. Saddle 
G(MN-—1)>0. 
b. Node 
Both G(MN—1) <0 
and 


1 
e [ary —— 4 At 0. 


Stable if S(G+1) >0. 
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c. Focus 


dod. 
G [acn - mea) ty 


Stable if S(G+1) >0. 


These stability conditions are summarized in Figure 3. 

It is evident from the stability diagrams that a variety of solutions is 
possible. Examples are given in Figures 4 and 5 which show solution curves 
for particular equations plotted directly with an analog computer. Coeffi- 
cients were chosen so that the isoclines are two parabolas intersecting in 
four points, one in each quadrant. In a biological application where vari- 
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Ficure 3. Stability of solutions near a singularity of the fourth group 


ables x and y represent populations, negative values probably have no 
meaning. A simpler set of curves results, however, with each singularity 
placed in a separate quadrant, and with negative values of x and y al- 
lowed. Proper choice of coefficients would permit a singularity in the first 
quadrant to have solutions near it similar to any of those of the examples. 

The only difference between the equations for Figures 4 and 5 is in the 
sign of coefficient a,, so that their isoclines and singularities are identical. 
The various solution curves result from different initial conditions, and 
time increases along the curves in the direction of the arrowheads. Nu- 
merical values of the quantities MW and G for each Group 4 singularity 
are plotted in Figure 6.°This figure is the diagram of Figure 3, and the 
kinds of solutions found in the examples are predicted accurately from it. 
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Ficure 4. Solution curves for « and y as implicit functions of time for a particular example, 
Time increases in the direction of the arrowheads. Different solution curves result from dif- 
ferent initial conditions. Seven singularities occur, the type of each being indicated. 
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Ficure 5. Solution curves for an example differing from that of Figure 4 only in the reversal 
of the sign of coefficient az. Singularities are located at the same points but the types are dif- 
ferent. 
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The general shapes of solution curves can be sketched from a knowledge 
of the location and types of singularities, of the isocline curves for y = 0 
and « = 0, and of the asymptotes for solution curves near a node or a 
saddle. Near either of these singularities, as time approaches plus infinity, 
the slope of a solution curve approaches the value (Farley, 1952) 

m= Ag. (16) 


Similarly, as time approaches minus infinity, the slope approaches the 


value 
A deo (eye) 


by = SS 


In these equations A, B, C, and D are coefficients in equation (10), A; is 
the more positive characteristic root and , is the less positive root, all 


6 G@=R/s 


e 
WI focus 
-2 


ineiablent SOG 0 stable 1f S>O 
' 


as 6. Stability diagram of Figure 3 with points located from examples of Figures 4 
and 5. Roman numerals identify each point. For points V and VII, S' > 0: f i 
aoe 5 ; for points VI and 


evaluated at the singularity in question. The algebraic sign of the slope of 
a solution curve reverses when it crosses an isocline for either y=Oor 
¢ = 0. A solution curve cannot leave the quadrant in which it starts. 
Using all this information allows solution curves to be sketched with some 
accuracy. 
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One type of solution of interest is that which is oscillatory, as near the 
focal points of Figure 5. The simplest pair of equations that give oscilla- 
tions are those studied by Volterra, 


: Az . 


y 


2 a 
Y=z (hy — Byx) y, 
y 


where a, is a negative quantity, but all other coefficients are positive. 
These are the same as equations (5) and (6), but with terms in x and y? 
omitted, and only the first terms of equations (7) and (8) used. It is not 
difficult to show that periodic undamped oscillations of small amplitude 
can occur near the one singular point (~, = ky/By, ye = Rz/B2), with 
angular frequency w = (—a,a,)”. 

If the terms in x” and +? of equations (5) and (6) are introduced in equa- 
tions (1) and (2), the effect is to produce damping. Damping, in general, 
modifies the frequency of oscillation and makes the solution nonperiodic. 
The amplitude either grows or decays, dependent upon the algebraic sign 
of the damping term. If terms of higher degree in the coupling functions 
of equations (3) and (4) are introduced into equations (1) and (2), the 
location of the singularity is changed and the frequency of oscillation is 
modified. Moreover, additional singularities may appear, and these may 
or may not have oscillatory solutions. 

In some oscillatory systems having nonlinear damping a limit cycle can 
occur. This is a periodic solution having its amplitude and frequency de- 
termined by parameters of the system, entirely independent of initial con- 
ditions. Small changes in the parameters do not destroy the periodicity. 
There seems to be no sure way of predicting whether a given system may 
have a limit cycle. It does not appear, however, that a limit cycle occurs 
with equations (5) and (6). In theory the coefficients might be adjusted so 
that the net damping is zero and a periodic solution would exist, as with 
the Volterra equations. The amplitude would then depend upon initial 
conditions and any small change in coefficients would destroy the periodic- 
ity. None of the examples studied has ever produced a true limit cycle. 

The work described here was done in part under research contract 
Nonr-433(00) between the Office of Naval Research and Yale University 


(Cunningham, 1954). 


The author is indebted to Professor G. E. Hutchinson for bringing the 
problem to his attention. 
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In a previous paper (Bull. Math. Biophysics, 16, 317-48, 1954) a transformation T of one 
graph into another was suggested, which may describe the relations between organisms of dif- 
ferent complexity. In this paper some topological properties of the transformation T are stud- 
ied. It is shown that the fundamental group of the transformed graph is homomorph to the 
fundamental group of the original graph. An expression is derived for the number of points in a 
point base of the transformed graph in terms of the number of points of the point base of the 
original when the point base of the latter consists only of residual points, and it is shown that 
the ratio of the number of points of the point base to the total number of points of the graph is 
in that case greater in the transformed graph than in the original. A combinatorial problem 
arising in connection with the transformation T is solved by deriving the number of possible 
ways in which » — n; indistinguishable elements may be arranged in 7; classes, permitting 
some of the 7; classes to be empty. 

The possible biological meaning of the increased ratio of the number of points of the point 
base to the total number of points of the graph is discussed. It is suggested that it may be 
interpreted as a decrease of regenerating ability with increase of differentiation of the organism. 
Those considerations suggest the possibility of deriving some general biological laws from the 
consideration of the properties of the transformation only, regardless of the choice of the pri- 
mordial graph. 


In a previous paper (Rashevsky, 1954) we outlined a topological ap- 
proach to general biology and discussed a transformation designed to 
describe some relational properties of the organic world. Familiarity with 
that paper is a prerequisite to the understanding of the following. For 
convenience we restate here the rules which specify the transformation T 
when biological functions are specialized in m classes of cells. 

Let n biological functions f, all arbitrarily numbered from 1 to n, be- 
come specialized in m groups of m, %%, . . . %m points each. We have Zn, =. 
n. Let group 7 specialize in ; biological functions f;, where / designates one 
of the m; points, and lose n — ; biological functions f,, where k designates 
one of the ~ — m; points. The m — n; biological functions f;, are divided 
into m — 1 groups, each consisting of m4, 2, .. . Mia, Miz1, .. . %m points 
respectively. We now formulate the following transformation rules 7: . 
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T,. Draw m distinct primordial graphs. The number m is a parameter of 
the transformation. 

T,. Remove from a graph i (i = 1, 2,...m) those m — n; points which 
correspond to its lost biological functions. 

T3. Connect each specialized point f\” of a graph 7 to those remaining 
points f of each graph & (k = 1,2,...i—1,7-+1,...m) which 
correspond to those points f\” of graph i to which f;” is connected, 
preserving the directions of connections. 


T,. Choose a distribution of m — 2; elements into 7; sets, so that the /th 
set has a number p\") of elements with 


Do) =n =H, (1) 


t= 


Some of the v‘")’s may be zeros. This will always happen when 2; > 
n — n;. This arbitrary chosen distribution plays the role of another 
parameter of the transformation. A third parameter is the choice of 
the n,’s. 

Ts. Inagraphz add — nu; subsidiary points, and divide them into 7; sets 
containing each v}’ points. Connect each of the v$” points of a set to 
all the other points to which the specialized point f\” is connected in 
the same way as f\”) is connected to them, preserving the directions of 
connections and also make the connections 


URAL isfaee ak UO: (2) 


Ts. Connect each point fi, (7 = 1, 2,...v{°) to a residual graph in 
the same manner as f\" is connected to its own residual graph. 

We shall now prove some theorems relating to the above transforma- 
tion. The primordial graph shall be denoted by P, the transformed graph 
by T(P). When reference is made to one of the m primordial graphs of 
which T(P) is composed according to T;, we shall speak of the component 
primordial graph to distinguish it from the primordial P. 

I. First we prove Theorem 1: The fundamental group of the transformed 
graph T(P) is homomorph to the fundamental group of the primordial graph 
Pi 

Theorem 1 may be considered as a direct consequence of the general 
theorem that in a continuous mapping of a complex onto another the 
fundamental group undergoes a homomorphic mapping (Seifert and 
Threlfall, 1934, p. 155). In the mapping of T(P) on P (loc. cit., pp. 334-35) 
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to each point of the graph T(P) corresponds uniquely a point of 7, and 
to each line of T(P) corresponds uniquely a line of P. This in itself does 
not guarantee a continuous mapping of T(P) on P, for the mapping of a 
line of T(P) onto the corresponding line of P may not be continuous. Li 
however, we assume, as we shall do, that it is continuous, then Theorem 1 
follows immediately. The following analysis, however, does bring out the 
more detailed nature of the homomorphism in question, which may be of 
some subsequent use biologically. 

As we have seen in loc. cit. we may visualize T(P) as consisting of m 
component primordial graphs, arranged in m layers, so that the corre- 
sponding points and lines lie above each other, and the corresponding 


FicureE 1 


specialized points, which we shall denote here by greek characters, a, 6, 
7 ..., are fused into one point. These fused points may be “pulled into” 
one of the m component primordials, for example, the top one. Each of 
these points a, 8, y . . . gives rise to connections (2) and is thus connected 
to v,, ¥g, etc., additional residual graphs.* The total number of the addi- 
tional residual graphs is 2(m — 1). ; 

Each point or line of a component primordial is mapped in a unique 
way onto a corresponding point or line of the primordial P. The same 
holds for each point or line of any additional residual graph. 

We define a path as a sequence of adjacent lines of a graph, regardless of 
whether the two adjacent lines have the same direction or not. A point 
following a path will follow, in general, some of its lines in the direction of 
that line, some in the opposite direction. Thus in Figure 1 the sequence 
01234 is a path; also the sequence 25067. 


* In our present notation connection (2) of loc. cit. would be written 


17 Ag. oe Ay, o 
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A path in which all lines are directed in the same manner is called a way. 
Thus in Figure 1 the sequence of lines 345681 is a way leading from the 
point 3 to the point 1. 

We now introduce the following notations. We denote by 


Wi (s,)5 Gim=ol ele ee tee An rr (3) 


a closed path that lies completely in the ith component primordial, origi- 
nates at the point p, and is characterized by a given sequence 5, of lines. 
In other words, to different s,’s correspond different paths. Thus in Figure 
1 the path 0123450 corresponds to s;, the path 05680 to sz, the path 
067810 to ss, etc. 

To each closed path W{?(s,) there corresponds a unique path W,(s,) 
in P as well as a unique path W{"(s,) in the kth component primordial. 

We denote by 


WEI (s); = Fs2,...,e(m—1);  r=1,2,...,@ (4) 


a closed path that lies completely in the z’th additional residual graph, 
originates at the point #, and is characterized by s, which may be different 
from the s, in (3). To the path (4) there corresponds in a unique way the 
path W,(s,) on P as well as the path W{"(s,) in the &’th additional resid- 
ual graph and the path W<"(s,) in the ith component primordial. 
We denote by 
Wao, (5) 


the open path that leads from the point a, to the point a along the line 
(2), connecting the two corresponding points a, and as. 
We denote by 


WE (s,) (6) 


the open path that lies completely within the ith component primordial 
and connects the point # to the point q, following a set of lines specified by 
s,. The path (6) corresponds in a unique manner to the path W,,,(s,) on P, 
as well as the path W{?(s,) on the Ath component primordial. 

We denote by 


ey (7) 


an open path that lies completely within the z’th additional residual graph 
and connects the points p and q, following a set of lines s,. To (7) corre- 
sponds in a unique way the path W,,(s,) on P and the path W“”(s,) on 
the k'th additional residual graph. 


SOME THEOREMS IN TOPOLOGY 115 


For convenience we shall use numbers for the subscript r in s,, thus 
writing Si, s2, s3...etc. Different indices correspond to different se- 
quences of lines. 

If a closed path originates at the point p, follows a set of lines s, to q, 
from there a set of lines s2 to r, and from r a set of lines s3 to p, then the 
closed path formed by the three open paths W%), W&, W{ is denoted by 
WS? (such isa ts Sa): 

We denote the fundamental group of T(P) by Gr, the fundamental 
group of P by Gp. The elements of Gr are classes of combinatorially 
homotop (Seifert and Threlfall, Joc. cit., p. 159) closed paths, originating at 
a fixed point. The choice of the point is immaterial in the present case 
because a change in this point merely results in an inner automorphism of 
Gr (Seifert and Threlfall, Joc. cit.). We shall therefore choose as such a 
point one of the specialized points, say a. The elements of Gp are classes of 
combinatorially homotop closed paths, also taken to originate in a. 

Now we shall determine the class of elements of Gr which map onto the 
identity element on Gp. 

First of all the identity element of Gr maps onto the identity element’ 
of Gp; for the identity element of Gr is represented by either the point a, 
which maps on the point a of P, or by any path which runs in T(P) from 
the point a to any other point (including a) and then returns to a retracing 
its steps. But to such a path there corresponds a unique path in P that 
starts at a, runs to any point (including a), and returns to a retracing its 
steps. But that path belongs to the identity element in P. 

We shall call a path that is combinatorially zero homotop an identity 
path and denote it by e. The identity element of the fundamental group 
we shall denote by £. 

To every closed path Wz in T(P) there corresponds a closed path Wp 
in P. Any closed path in T(P), combinatorially homotop to Wz, is ob- 
tained by a) either following or preceding W,r by identity paths, or b) by 
extending Wr at some point with a closed path W* which starts at that 
point pr, runs to some point gr, and then is retraced back to pr. Since to 
any identity path in T7(P) there corresponds an identity path in P, there- 
fore to any path in T(P) which is combinatorially homotop to Wz because 
of a) there corresponds a unique path in P combinatorially homotop to 
W p. Since to any path W7 there corresponds a path W>, which starts at 
some point fp, runs to a point gp, and is retraced to pp, therefore to each 
closed path in T(P), which is combinatorially homotop to Wr because of 
b), there corresponds a path in P which is combinatorially homotop to 


Wp. 
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Hence if one representative path of an element of Gr corresponds to a 
representative path of an element of Gp, then to every representative path 
of that element of Gr there corresponds a representative path of the ele- 
ment of Gp. 

The path in T(P), 


W® (s,) [WS CS) as t,kR=1, 25. -5 5% (8) 


does not belong to the class which is the identity element in Gr unless 
i = k. But since WS (s,) and W(s;) both map onto W,(s:) in P; there- 
fore, (8) maps on a path of the identity element in Gp. 

By a similar consideration we see that 


W,,.WE (5,) Wey WE (5) 1 Wg) 11,1 (9) 


b 

which does not belong to the class which is the identity element in Gr 

unless 7’ = k’ maps onto a path of the identity element of Gp. (Remember 

that a, and a» as well as the paths Waa, and Wa,«, all map on a in P.) 
The path in T(P) 


WO (5) W,, We? (s,)17We (10) 


for the same reasons maps onto the identity path in P. 
Another non-identity element in Gr which maps onto the identity ele- 
ment in Gp is represented by 


WO CS LW Sa) la (11) 


because the point moving along this path maps on a point which in P 
moves along W4,(s1) Wg (51). 

If p is a point of P adjacent to a specialized point a, p; the corresponding 
point in the zth component primordial graph, and Wa,», the line joining 
a, to p; according to T;, then the element in Gr, represented by 


UES: [Weel bake (12) 
also maps onto the identity element in Gp, because Wa,«, maps onto the 
point a, while both Wa», and Wa», map onto Wap in P. 

All elements of Gr represented by paths of the form (8), (9), (10), (11), 
(12), and their products map onto the identity element of Gp. 

From those elements we can obtain others that also map onto the iden- 
tity element in Gp, by “intercalating” other paths of the form (8) to (12). 
For example, between the two paths in (8) we may intercalate a path 
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W (s2)[W(s2)|-!, which is also of the form (8). We obtain 
ee OSS) IW) Cs, Weis) =e (13) 


to which in P corresponds 


Vaeese) Wes») Ww, (s,) | =a Wes) ] =] — e. (14) 


We shall call the element represented by the path (13) an extension of the 
one represented by (8). The element represented by (13) may again be 
extended, etc. 

We may extend an element represented by a path of any form (8), (9), 
(10), and (11) by intercalating other paths of any of these forms. For ex- 
ample, from (11) we obtain, by intercalating a path of the form (10), 


Wee (8) WP Cs.) Weg, (WD? (5,) 1-1 Woes (WO (s,)] 2, (15) 


which, since Ws, maps on the point 6 in P, maps in P onto 
Wp (5) Wal(syW,* (s,) Wy s,) =W (5) Ws) =e. (16) 


Allelements represented by paths of the form (8) to (12) as well as their 
products and extensions exhaust the class of elements of Gr that map onto 
the identity element of Gp. This is seen from the following considerations. 

The paths of the form (8), (9), and (10) are all constructed in such a 
way that two corresponding closed paths are followed in opposite direc- 
tions. Those two corresponding paths may lie in the following three man- 
ners: 

1. One in a component primordial graph, the other in another com- 
ponent primordial graph [(8)]. 

2. One in an additional residual graph, the other in another additional 
residual graph [(9)]. 

3. One in a component primordial graph, the other in an additional 
residual graph [(10)]. 

The paths of the form (11) represent closed paths that lie in two dif- 
ferent component primordial graphs in such a way that the two open 
paths, of which they are composed and which connect two specialized 
points, are corresponding paths followed in opposite directions. 

Finally there remain the paths of the form (12) that are due to connec- 
tions Ts and ioe 

One might have thought of the possibility of having one closed path in 
a component primordial graph, and the corresponding closed path, fol- 
lowed in the opposite direction, in an additional residual graph that be- 
longs to subsidiary points other than those of a. This is impossible because 
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the other closed path in the additional residual graph must be reached in 
that additional residual graph from a point which corresponds to a. But 
only the additional residual graphs that belong to the subsidiary points of 
a have such points. 

Now consider any other element in Gr that does not map onto the 
identity element of Gp. Take, for example, an element S represented by a 
path of the form 

Wo (s,) WH (5,) W,, WEY (5) WWD (5,) Wi (s,) WO (s,). (17) 


a 


This element maps in a unique manner onto the element S* represented by 
Wa (51) Wa (52) Wa (53) Wa (Se Se Se) (18) 


By giving toi, k, p, g, andr in (17) all possible values from 1 to m, to 7’ 
all possible values from 1 to w(m — 1), and my multiplying them by ele- 
ments represented by paths of the form (8), (9), (10), (11), and (12) or by 
any products of those elements we find all elements S in Gr that corre- 
spond to the element S* in Gp. 

Take another element Si, represented by 


Wi s,) We WEN SY WOKS,) Ws oe 


To it there corresponds the element Sj‘, represented by 
Wa (57) Wa (53) Wa (59). (20) 


The element SS; of Gr is represented by the path in T(P) followed by first 
traveling the closed path (17) and then the closed path (19). The represent- 
ative point in P will first travel the path (18) and then the path (20). 
Hence if S corresponds to S* and S; to S¥, then the product S'S; in Gr cor- 
responds to the product S*S¥ in Gp. Also to S~! in Gr there corresponds 
[S*]-* in Gp. Thus we have established a homomorphic mapping of Gr 
onto Gp, with the elements represented by paths of the form (8), (9), (10), 
(11), (12), and their products as the kernel of the homomorphism. 

It now follows from a general theorem in the theory of groups (Pont- 
rjagin, 1946, p. 11, theorem 1) that all elements of the form (8) to (12) and 
their products form a normal subgroup of Gr and that all elements S that 
map on the same S* are cosets of that normal subgroup. 

As an example consider in Gr the element S’, represented by 


W (s,) WO (s,) WO (s,) (21) 


which maps onto the element S’*, represented by 


, Wa (51) Wa (52) Wa (83) (22) 
in Gp. 
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Now consider the element S** in Gp, represented by 
[WE (s,) 17 TW (s ) 1 [WO (s,) 17 WO (s,) WO (s,) Wi) (s,). (23) 


It maps onto the identity element E in Gp, being obtained from elements 
represented by paths of form (8) by two extensions. Multiplication of .S’ 
by S** gives an element S’S** of the left coset, represented by 


W®) (s,) W® (s,) WO (s,), (2.4) 


which again maps onto the element S’* of Gp. Multiplication of S** by S’ 
gives an element S**S’ of the right coset, represented by 


[WO (s,) 1 W® (5,)]-1 [WO (s,)] 1 W® (s,) WO (s,) WO (s,) 


| (25) 
XW (s,) W® (s,) WO (s,), 


which also maps onto S’* in Gp. 

In Joc. cit. we found that the transformation T is too general from a 
biological point of view, giving many more different possible organisms 
than actually occur. We suggested on page 342 that the requirement of 
homomorphism of the fundamental groups of T(P) and P may perhaps 
impose a restriction upon T(P). As we see now, the fundamental group of 
T(P) is always homomorph to that of P for any choice of the parameters 
in J. Hence other restrictive conditions must be looked for if we are to 
retain this particular transformation T at all. As remarked on page 325 of 
loc. cit., it is not likely that this particular transformation will be of bio- 
logical use. However, since transformations of this type have not been 
studied hitherto in topology, to the best of our knowledge, it becomes 
necessary to first study im abstracto topological properties of all sorts of 
conceivable transformations if we are to eventually find the most useful 
one. 

II. We shall now study some relations between the point bases of P and 
T(P). We remind the reader of the definition of a point base (K6nig, 1936, 
p. 88): 

A set B of points of a graph K is a point base of K if 

1. To each point # of the graph, not belonging to B, there exists a way 
from a point of B. 

2. There is no way from any point of B to any other point of B. 

The point base Bp of P may consist either of residual points only or of a 
specializable point only, or some points of the point base may be residual, 
some specializable. We shall discuss here only the first case. 
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Lemma: If in a graph K we remove a point p, which does not belong to 
the point base B of K, then the point base B’ of the partial graph K’ thus 
obtained consists either of the set B or of the set B and of some additional 
points. 

Proof. If any point of the graph K can be reached from a point of B by a 
way that does not pass through # then the removal of p does not alter the 
point-base properties of B and B remains the point base of K’. Suppose, 
however, that there are w points that can be reached from a point of B 
only by ways that pass through p. Denote by A, the set of these w points 
and by K,, the partial graph of K that consists of these w points and of 
their connecting lines. Since K,, is a directed graph it has a point base B,, 
(Konig, loc. cit., p. 89). We assert that the point base B’ of K’ is the totality 
of the points belonging to B and B,, which we shall denote by B + By. 
The first point base criterion is satisfied because every point of K’ can be 
reached by a directed way from a point of B + B,. Any point of K’ not 
belonging to A, can be reached by a directed way from a point of B, 
while any point of A, can be reached by a direct way from a point of By, 
by definition of B and B,. The second criterion of a point base is also 
satished by B + B, because by definition of B and of B,, there exists no 
directed way in K’ which connects any two points of B or any two points 
of B, and, because of the removal of the point p from K, no directed way 
connects any point of B with any point of B,,. Hence the point base of K’ 
is either B or B + B,, which proves our lemma. 

We introduce the following notations: we denote by p the point in the 
ith component primordial of T(P) that corresponds to the point p of P and 
by p® the point in the 7’th additional residual graph of T(P) that cor- 
responds to p in P. By Ra, we denote the additional residual graph that 
belongs to the subsidiary point a, of the specialized point a. If A isa set of 
points of P, we denote by A“ and by A“ the sets of corresponding points 
in the ith component primordial graph of T(P) and in the 7’th additional 
residual graph of T(P) respectively. We now prove 

Theorem 2: If Pp is a point base of P that consists of residual points 
only, then a point base of T(P) consists of all points corresponding to the 
points of Bp in every component primordial and in every additional 
residual graph plus such points in each Ra, that correspond to the points 
that are added to the point base Bp of P in order to form the point base of 
the partial graph of P, obtained by removing from P all specializable 
points except a. 

Denoting by C the set of all these additional points and by Br the point 
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base of T(P), we may state the above in symbols thus: 


Era erat > BOL, (26) 


the summations being extended over all z’s and i” ’s. From the lemma it 
follows that C may be the null set, though, in general, it is not. 

Proof. a) Any point of a component primordial graph i in LP ican 
be reached by a directed way from a point of B$, since the ith component 
primordial is identical to P, having all specialized points, because of 7. 

b) An additional residual graph Ra, consists of all residual points plus 
the point a, that corresponds to a in P. Hence its point base consists of 
BE” plus these additional points that must be added when all specializable 
points except a are removed from P. These additional points are part of 
the set C. Every point of Ra, can be reached from a point of that base by a 
way. 

c) Any point a,, B;, etc. can be reached from a point of any BS by a 
directed way by first reaching in the ith component primordial the point 
a by a way that is possible because of a) and then following the directed 
way aa,, etc. 

Thus the first criterion of a point base is satisfied for the expression (26). 

No point of BS” can be reached by a way from any point of BY or of 
B&”. The only manner to reach by such a way a point of BS” from a point 
of BY is to follow in the kth component primordial the directed way from 
some point of B£” to one of the specialized points, say a, and then follow 
a directed way in the ith component primordial from a toa point of Bf”. 
To such a way there would correspond a directed way in P from one point 
of Bp to another, which is contrary to the point base properties of Bp. 

In a similar manner we see that there exists no way from any point of 
BY? to any point of BY. 

It remains to prove that there exists no way from any point of B® to 
any of the points contained in the set C, nor from any point of B ° to any 
points of the set C. 

If there were a way from a point p© of BS to a point q of the set C, 
this way would have to pass through a specialized point a, then follow the 
way aa,, and then from a, reach the point g@. Since aa, maps on the point 
a in P, therefore, such a way would map on a way in P that goes from the 
point p of Bp to the point g outside of Bp, passes through a but does not 
pass through any other specialized point, because in the z'th additional 
residual graph R., there are no points corresponding to those other spe- 
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cialized points. This is, however, impossible because the point g is one of 
those which can be reached in P from Bp only by way of one of the re- 
moved specialized points 8, y,.... It is their removal that added the 
point g to the point base of the residual graph Ra. Thus the theorem is 
demonstrated. 

Definition. The ratio of the number of points of a point base to the total 
number of points of a graph is called the point base ratio. 

Corollary of Theorem 2: The point base ratio of T(P) is greater than 
the point base ratio of P. 

Proof. Let the number of points of P be a) and number of points of Bp 
be po. Then the point base ratio of P is 


Yaseen (27) 


On each residual point of P there are mapped m + n(m — 1) points in 
T(P) (loc. cit., p. 335). If p, denotes the number of points in the set C, then 
the number of points in Bz is po[m + n(m — 1)] + p. while the total num- 
ber of points in T(P) is aj = may + n(m — 1)(ao — n) [loc. cit., eq. (5)]. 
Hence the point base ratio rz of T(P) is 


_ polmtn(m—1))] the _ polm+n(m—1)] + 2. 
maytn(m—1) (a—n) aolm+n(m—1)] —n?(m—1) 
(28) 
gsi 
Qo 


the inequality holding even for p, = 0, because m > 1,” > 2. 

III. Two combinatorial problems are connected with the transforma- 
tion T: one is the question of the total number of possible ways of dis- 
tributing the ~ distinguishable specialized points among the m component 
primordial graphs so that each of the m component primordial graphs re- 
ceives at least one point. Biologically this means the question of the total 
number of possible ways of distributing » specialized biological functions 
among m tissues. The expression for this number has been derived else- 
where (Rashevsky, 1955). 

The other problem that we shall now treat is one of the possible number 
of distributing ~ — m, non-distinguishable elements into ; sets, empty 
sets being permitted. This problem is posed in connection with T,. 

The problem is somewhat similar to that of partitioning an integer p 
into the sum of q < p integers except that in the latter case all g numbers 
must be non-zero. The number NV) of all possible partitions of ping 
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classes is (Netto, 1901, p. 11) 


Se ee 
OE NGe i): (29) 
We shall denote by Q® the number of possible ways of distributing m — n,; 
elements into 7; classes, classes with no elements being permitted, and we 
shall prove that 


n—1 


Opes emis) (30) 


Equation (30) is correct for any » and m; = 2. In fact, all possible ar- 
rangements of ~ — 2 elements into two classes are exhausted by the fol- 
lowing scheme 


n—2, 0 
n—3, 1 
n—4, 2 
b 
Ome n—2 


Altogether there are 2 — 1 = (3-}) possibilities. We now shall prove 
that if (30) holds for any and for a given ;, it holds also for the same 
n and for m; + 1. 

We now have to distribute ~ — 1; — 1 elements into 1; + 1 classes. 
All possibilities can be exhausted in the following manner: 

We first distribute the  — n; — 1 = (m — 1) — n; elements in all pos- 
sible ways into ; classes, leaving the (7; + 1)st class empty. The total 
number of such possibilities is 


(n—1) 
OT - 


Next we put one element into the (#; + 1)st class and distribute in all 
possible ways the remaining (m — 2) — m; elements in the m; classes. The 
number of such possibilities is 


0629 


Thus we proceed until we have m — n;-—1= (n — 1) — n,elements in the 
(nm; + 1)st class and none in the first ; classes. The number of such possi- 
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bilities is 
Qe =Qu=1, 


according to (30), since (30) holds for m;, and as is also directly clear. 
Hence we have 


a=N—ni 


Os vii >» rea (31) 


Since (30) holds for 2; and any ”, therefore, substituting (30) into (31) and 
using a well-known identity [Netto, 1901, eq. (10)], we find: 


Onna (n= 1) ean) om oils 
Q.E.D. 


IV. The corollary from’Theorem 2 suggests an interesting possible bio- 
logical interpretation. As we have pointed out in Joc. cit., p. 327, a biologi- 
cal function f represented by a point of the graph does not take place, 
unless the “‘preceding’’ biological functions, the representative points of 
which are connected to f by arrows directed toward f, also take place. 
Hence if we remove or inhibit all biological functions, the representative 
points of which constitute the point base of the graph, the organism will 
not be able to function at all. 

Now let us make an additional hypothesis, namely, that an impaired or 
removed biological function will “regenerate”’ if all its “preceding” bio- 
logical functions are intact. Then an impairment or removal of some bio- 
logical functions of an organism by impairment or removal of some ana- 
tomical parts of it will result in a regeneration of the biological function, 
and, therefore, of the anatomical part which carries it, as long as the 
biological functions which correspond to the points of the point base are 
intact. 

If the point base ratio is very small, then the chances of impairing or 
removing the biological functions corresponding to the point base, by ran- 
dom removal of other functions, are small. The probability of regeneration 
of lost biological functions is large. On the other hand, when the point base 
ratio is large, and near 1, then the probability that a random removal of 
any biological functions also removes some of those which correspond to 
the point base is great, and the probability of regeneration is small. 

With the above interpretation, the corollary, translated into biologi- 
cal terms, thus states that if the point base of the primordial graph con- 
sists of only residual points the more complex and highly developed and dif- 
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ferentiated organisms show a lesser chance for regeneration of biological func- 
tions lost by injury than do the lower organisms. This is a well known general 
biological fact. 

Thus stated the theorem also emphasizes that we deal only with prob- 
abilities of regeneration. The probabilities of regeneration of an organ or 
biological function in the same organism will vary from one biological 
function to another, depending on the exact topology of the graph. But the 
general statement italicized holds regardless of the particular choice of the 
primordial graph, as long as the corollary holds. We do not wish to put 
too much weight on the above interpretation of the corollary, this inter- 
pretation being susceptible to several criticisms, but we believe that the 
following important possibility emerges from the above discussion: It may 
be possible to derive some general biological laws by considering only the 
properties of the transformation, regardless of the choice of the primordial 
graph. We may then formulate the following problem: What conditions 
should a transformation satisfy in order to lead to given general biological 
laws? 

Is there any reason to believe that biological functions which correspond 
to the point base of the primordial graph are all of “residual” nature? It 
seems that such biological functions are actually found in the biological 
functions of activities of the genes. As we remarked in loc. cit., the pri- 
mordial graph discussed there is undoubtedly highly oversimplified, and 
the actual graph is much more complex. The hundreds of biological reac- 
tions which take place in a cell are all in some way determined by the 
presence of appropriate genes. Thus, without actually knowing the physio- 
logical nature of the action of the genes, we may readily identify the rep- 
resentative points of those actions with the point base of the graph. The 
equal distribution of genes upon the daughter cells in mitosis expresses the 
fact that the representative points are residual points. Possible interac- 
tions of some genes may seem to violate the second condition of the point 
base, indicating the existence of directed ways between some representa- 
tive points, which correspond to genic activities. In such a case, however, 
though the set of representative points of the activities of all the genes do 
not form a point base, it is readily seen that this set of points has its own 
point base which is also the point base of the primordial graph, and is still 
composed of residual points only. 

In the primordial graph shown in Figure 1 of Joc. cit. the point base ac- 
tually consists of specializable points. As remarked, however, that graph 
is used only as a tentative illustration. An appropriate biologically plau- 
sible change in it easily avoids the difficulty. 
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The gene functions, according to this picture which is to be considered 
primarily as an illustration of a possible line of thought, correspond to 
the point base of the primordial graph. According to the Lemma, any 
specialization of tissues, connected with loss of some biological functions, 
results in general in the circumstance that points which correspond to 
other than gene functions are added to the point base of the organism. 

We must, of course, keep in mind that the genetic constitution of an 
ovum, from which a higher organism develops, is different from the genetic 
constitution of the primordial cell, whichever this turns out to be. The 
former genetic constitution is much more complex. It also varies from 
organism to organism. This, however, does not necessarily invalidate the 
above interpretation. The genetic constitution of the primordial organism 
may change, the number of different genes increasing. Only when enough 
“new” genes are available to cause the differentiation of some function 
does a higher organism develop. The relations required by our interpreta- 
tion of the corollary still hold in this case. 


The author is indebted to Dr. Ernesto Trucco for a thorough discussion 
of this paper. 

This work was aided by a grant from the Dr. Wallace C. and Clara A. 
Abbott Memorial Fund of the University of Chicago. 
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Simple linear systems in which there is a steady-state velocity may under some conditions 
have solutions which are periodic but damped. However, even under favorable circumstances 
it is not likely that the periodic property can be differentiated from a single overshoot because 
the effect of damping is so great. 


It has been shown by J. Z. Hearon (1953) that linear systems, in which 
the steady state is also the equilibrium state, cannot show periodic solu- 
tion. If the system is, for example, an enzyme system in which there is a 
steady-state velocity, the conclusions no longer hold. It is the purpose of 
this note to investigate some aspects of this problem, since in most bio- 
logical systems one is concerned with systems having steady-state 
velocities. 

Consider the following system of reactions 


hy 
M+ SH (X4S1) 


ke 


(X19) a Xetli 


- (1) 
DO So (X2S2) 


If we may ignore the complexes this reduces to the following, if there 
are only 3 X’s, 
XitS1= XitPi, 


XotSo x XitPrs, (2) 


Xot+Ss = XitP3. 
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On the other hand, if the S; and P; are large enough the first set of reac- 
tions reduces to 


(X Soe = (X22) , 


(XoSe) = = (XsSs) , P (3) 


(X35 = a 51) - 
If in (2) S; and the P; are present in sufficient quantities so that over the 
time considered they do not change appreciably or are maintained nearly 
constant, then if k; = mS and k_; = m_,P;in (2) and k; = n;in (3), etc., 
we find that the linear differential equations for both cases are as follows 


b= — (RRL, 4) eRe TR Bias» (4) 


where «; = (X,) in (2) or x; = (X,S,) in (3) and the subscripts are mod 3. 
The characteristic roots of (4) are 


N= 0 
1 ! ! rape 
= ayes (ki+k!) [14 Vi-w, (5) 
where 


A(Rik3+RiRs+RSRi +R Rk ot hiok sth gk ytkikiothkgk_3t+k3k_y) 


JS ae yt a 


It can be shown that w cannot exceed $ for non-negative values of k4,. 
Now y has the maximum value 4 when k = k, = ki = k’ andk!, = Oor 
when k; = 0 and k!, = k!, = k!, = k’. These correspond to the cases in 
which the rates in one direction in the cycle are negligible compared to 
those in the reverse, the latter being about equal to each other. This case 
then gives imaginary roots and the frequency relative to the rate constant 
22(k1.; + k1,) of the damping factor, ~/y — 1, has its greatest value. 
Since one is interested in the possibility of roots which can be detected, 
it is necessary to determine the maximum effect of the periodic term. This 


can be estimated from the amplitude at the end of one complete cycle. 
Thus when 


RVpy-1) t= (Vi) Rt=H27, 
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the amplitude is reduced by a factor 
€ Cai a= es = 00002 . 


From this it can be seen that although the solution has imaginary time 
constants, the periodic character of the time course of the changes could 
not be detected under ordinary circumstances. For the case in which 
A — B—C-—A and for the case in which initially B = C = 0, then A 
undershoots by about 23%, B overshoots by about 15%, and the over- 
shoot of C is negligible. The amount of subsequent over- or undershoot is 
negligible. Thus the most favorable case for periodic behavior would not 
ordinarily be distinguishable from a single overshoot. 

For the case of four steps the situation is more complicated. If we con- 
sider the situation corresponding to the maximum periodic effect for the 
three-step case we can show that the cyclic case in which the rates in one 
direction are negligible compared to those of the opposite direction, the 
latter being equal to each other, corresponds to a maximum in the follow- 
ing sense. The ratio of the imaginary part to the real part of the complex 
root decreases whenever one of the rates is made different from the others. 
On the other hand, introducing a non-negligible rate in any one of the re- 
verse reactions causes a decrease in the ratio. The introduction of a nega- 
tive rate causes an increase in the ratio but we consider here only that 
positive rates are physically meaningful. In this case the ratio is again so 
unfavorable that it is unlikely that the periodic effect can be differentiated 
from a simple overshoot. Hence, for the system considered, although peri- 
odic solutions are mathematically possible, we find that the effect of damp- 
ing is so great that the periodicity is not likedly to be observed. 

This work was aided by a grant from the Dr. Wallace C. and Clara A. 
Abbott Memorial Fund of the University of Chicago, as well as by re- 
search grant G-3312 from The National Cancer Institute, of the Nation- 
al Institutes of Health, Public Heath Service. 
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The effect of competition between species is considered in terms of a mathematical model. 
It is found that, except for special relations among the parameters, only one species should 
remain, as is found experimentally. It is also found that unless new factors arise in the competi- 
tion, one can calculate an index to predict which of the two species is most likely to become 
extinct. The index involves quantities measured from isolated populations. A preliminary 
extimate yields satisfactory values when compared with the experimental results of competi- 
tion. 


In a previous paper (Landahl, 1955) equations were derived to repre- 
sent the time course of the population numbers of various stages of the 
flour beetle. Thus expressions were obtained for the number of eggs (£), 
larvae (L), pupae (P), and adults (A) as functions of the time ?¢, the 
expressions involving survival functions for the various stages. The sur- 
vival functions in turn were supposed to depend upon mortality due to 
the interaction of the various stages. It is the purpose here to consider the 
effect of competition between species in terms of this model. 

Competition between two similar species, no differential resistance to 
cannibalism. Consider the situation in which two kinds of insects are in- 
troduced simultaneously into a container. We shall suppose that each kind 
has been studied separately under identical conditions. Suppose, further, 
that there is no interbreeding and that there is no differential resistance 
to cannibalism, that is, a given stage of either type cannibalizes another 
developmental stage to the same extent whether the latter type is the 
same or different from its own type (cf. Birch, Park, and Frank, 1951). 
Consider also that there is no interaction between any mortality factors. 
In this case one can write two sets of equations, one for each type, but in 
which the values of the survival functions are the same for both kinds. 
The logarithm of S; of (28) of Part I (Landahl, Joc. cit.) is now just the 
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sum over j of Cj; for one species plus the corresponding sum for the 
second. 

If, for example, both kinds have the same times for the duration of the 
various stages and the same death rate for adults, and if we start with 
adults only, with the same number of each kind, then for the duration of 
the first life cycle, the various stages will be in the proportion to the egg- 
laying rates of the two species, Ry:rg. During the next life cycle time the 
added adult populations will be in the proportion R%:rz if the initial 
number of adults is small compared to the steady-state value, and if we 
assume that the egg-laying rate is inherited exactly. After this the ratio 
between the numbers of the adults added per unit time changes in the 
same direction but more slowly. As a result of the death rate the total 
population becomes limited so that the increased ratio between the two 
species means a decline of the lesser with an increasing probability of 
extinction. A similar argument shows that if two species are different only 
with respect to death rate (D or d) or to the duration of one of the develop- 
mental stages, one of the species will become extinct. Differences in degree 
of cannibalism, when undirected, will not alone produce extinction of one 
of the species. 

When several factors are different it is easiest to see the result if we 
consider the situation in which two populations, which are each in a stable 
state, are mixed in a given proportion. If, for example, the proportion is 
one to one, each population can be divided in two at the census so that 
each mixture has the same volume as originally. We shall consider the 
possibility of a factor J which results in a constant death rate in the larval 
and pupal stages. Let S; represent the survival function for this effect 
acting over the entire period over which it acts. The survival functions for 
the egg stage for the small larval stage, etc. are denoted by S72, Ss, etc., 
the quantities Tz, T's, ... denoting the durations of the egg stage, the 
small larval stage, etc. (Landahl, Joc. cit.). Thus Sy, Ss, . . . represent 
survival probabilities for unit time. Let one population be represented by 
capital letters and the other by lower case letters. Then since AA and Aa 


are approximately zero in the separate populations, we have (Landahl, 
loc. cit.) 


RES PSPSSE*S USS = 5 7 stasisstx stu sizes = d , (1) 
In these expressions it is understood that the stable values of the respec- 


tive population numbers are used. If we mix volumes of the stable popu- 
lations in the proportion P to p (p + P = 1), we will find that AA and 
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Aa will no longer be zero but given by 


AA =P4aD/| x)". - =) 'P- 1 


Ra d SJ Dp 
=PA D\t —= — = 5 fates Ts—tg ps Tut Tr—tz > Tp—t = 
(oe) Yn D oF Si Sg Sur 55” Sp P 1 


(2) 


with Aa given by expression (2) with capital and lower case letters inter- 
changed, except in the subscripts. If the various values of the survival 
functions can be estimated from equations (21) through (30) of the pre- 
vious paper, then in the absence of differential interaction one can predict 
the rate of change for each population. Note that since generally the 
values of the various survival functions will be at least comparable stage 
for stage, then to this approximation the terms s7*~‘= and Sts"? = 
1/S*2“2 are roughly reciprocals. Furthermore, where T, is nearly equal 
to t,, the term in s, drops from (2). Unless some stage has rather different 
stage durations and also very different values for the survival coefficient, 
the term under the exponent # in the expression (2) for AA is approxi- 
mately the reciprocal of the term under the exponent P in the expression 
for Aa. Hence if one term is greater than unity, the other is likely to be 
less so that AA and Aa are likely to be of opposite sign, one population 
increasing, the other decreasing, but not necessarily at the same absolute 
rates (cf. Park, 1948, p. 303). 

From the discussion above it is clear that when the ratio between the 
terms which have the exponents P and p in (2) and in the corresponding 
expression for Aa is large, the population represented by the capitals is 
favored. Hence, if we form this ratio denoted by 3? we have an index which 
is perhaps the simplest single index from which to predict the outcome, 
values greater than unity favoring the population represented by capitals. 
The ratio, $, may be shown to be equivalent to 
Ss @ pra (se Piet Pte Aa £21)5,) aFe(S (eFs... 

Sz Sp re D sy (3) 


a = ‘R d ‘iS. — pal FT ere eee ATs 2(€4—€,) 
X Ga (5S (R ,) 1 (D) :(S,) 4 re (Ae 


where bar values are averages of the form S = s?S'-?, 4 = (Tx — tz)/Tx, 
etc., it being assumed that the e’s are sufficiently small so that e? can be 
neglected compared with e. The third expression involves the use of ex- 
pressions (4) and (5) below, and so holds only for the special case to be 
considered next in which all larvae are considered to be equivalent. Note 
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that the proportions of the species in the mixture enter only through the 
averages 5 and hence the results depend only slightly upon the proportion 
P. However, probability considerations indicate that the dependence on 
the ratio may not be neglible, especially when the numbers involved are 
small. 

As an example we shall consider the case of the interaction between 
Tribolium castaneum Herbst and T. confusum Duval studied extensively 
by Park (Joc. cit.). Since only the totals for larvae and pupae are given and 
since we shall here discuss quasi-steady states, we shall set S's8Sy"S7z? 
equal to S7+ where now T, is the sum of the previous subdivisions and 
S7: is now the survival function over all larval stages. Also we shall sup- 
pose that mortality not due to interaction is confined to the pupal stage. 
In this case we find for the steady state 


eLLEEA 
FE eae 4 
on SNe ie (4) 
ess 
Tr = Bs 5 
Sz R , A2,DT?,’ (>) 


from which we may estimate the mortality effect for eggs and larva from 
the steady-state values of each separate population. It should be kept in 
mind that the lumping of the mortality effects of all larval stages is apt 
to give rather poor values. However, for comparative purposes the result 
may be satisfactory. 

From the data given in Park’s paper we have the following values for 
T. castaneum (capitals) and T. confusum (lower case) when neither is 
parasitized: Rz, rz, D, d are respectively 6.4, 5.7, 0.0009, 0.005 per day; 
Tz, tz, Tx, tr, Tp, tp are respectively 3.8, 5.5, 22.8, 22.4, 6.2, 7.0 days. 
For the values of the population numbers in the steady state we use those 
corresponding to non-parasitized cultures: A,,, @., L~, l.. are respectively 
15.7, 12.7, 17.6, 15, 6.2 individuals per gram, assuming that the pupae 
and larval values represent mostly larva. We shall assume that the egg- 
laying rates and the death rates are those under steady-state conditions. 
We also assume that S; = sy = 1 in the absence of infection and neglect 
pupal cannibalism. From the above values we estimate that 


STt = 0.034 , st= 0.053, Si = 0.042 , s#= 0.0166. 
Noting that T;, = t,, we then find for the 
AA = 0.142P [2.2!-P — 1] = 0.009 [2.2!-P — 1] A 
Aa= — 0.064(1—P) [1—0.39?] = — 0.005 [1 —0.39?] a ; } ~ 
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P being the fraction of T. castanewm. Thus, regardless of the fraction as 
the T. castaneum adult population increases while the other decreases. If 
P= 1/2, AA = 0.035 and Aa = —0.012 = 0.00194 so that in the ab- 
sence of transient effects it would take about 100 days for the adult num- 
ber of T. castaneum to reach half way toward its stable value in the ab- 
sence of T’. confusum, if the rate remained constant; for T. confusum the 
time would be about 300 days to drop to one half. However from (6) it is 
evident that the former rate of increase would decrease to zero while the 
latter decrease would become more rapid on a relative basis, approaching 
—0.0031a. Since the data by Park (oc. cit.) show that T. confusum does 
become extinct in 12 out of 18 experiments, the above results are not in- 
compatible with the data. 

A similar calculation could be carried out for the interaction in the 
presence of the Adelina infection if the effects on the egg-laying rates, 
death rates, etc., were known. It may be pointed out, however, that from 
the values A, .../,., equal to 4.2, 9.1, 7.9, 9.9 respectively (averages of 
I, II, and III of Table 10, Park, 1948) one finds the quantity (A,,/L,,) 
in (3) is reduced by the factor 1.6 due to the Adelina infection. Hence even 
without any other bias against T. castaneum it is only necessary for S;/s; 
to be less than about 0.6 for T. confusum to be favored. Actually, adults 
may possibly die from the infection so that the pupal (or very late larval) 
death rate due to the infection does not even have to be as severe as 0.6. No 
estimates were found but this restriction seems compatible with the state- 
ment that deaths owing to this cause were plentiful. 

Having estimated the survival functions in the steady state we may at- 
tempt to reconstruct each population history from these values together 
with the given data, assuming only simple interaction. Since only a linear 
combination between A and L is given from S, in the stable state, it is 
necessary to make some assumption about the relative cannibalistic tend- 
encies of the adults and the larvae. Since 7. castanewm is more variable and 
the larval counts do not drop greatly when the adult population reaches 
its maximum, and in view of the data discussed in Figures 1-5 (Landahl, 
loc. cit.), we shall suppose that the larval pressure is not very much less 
than that of the adults. We choose the ratio of 1/2 somewhat arbitrarily 
for both egg and larval interaction. Then for 7. castaneum we find 


S72 = exp {—0.129(A4+3L) } 


(7) 
STi = exp {—0.137(A+3L)}. J 
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These expressions are then introduced into the expressions (25) which 
for this case reduce to 


E(i) =24A (¢— 0.4) VS2#(¢— 2) , (8) 


L@ =144A ¢—3.35) S@U@—2.9) 4/52 iia 2, (9) 
AA ($— 1). =29A G—8.39):ST8(¢— £9) S75 er 08 Se (10) 


For T. confusum we must choose a much larger ratio between the adult 
and larval cannibalism. One reason is that the populations are less variable 
so that in view of the previous discussions (Landahl, Joc. cit., cf. Fig. 6) 
the adult pressure is probably large. Also, since the larval populations de- 
crease to very small values when the adult population overshoots in the 
early history of the population, a much larger ratio is required. Again, 
since no data on egg numbers are available, the same ratio is used for both 
Sy and Sz, a value of 1/10 being chosen somewhat arbitrarily. Thus for 
T. confusum we have used 


sg=exp {—0.308(4+ 5 L)} } 


(11) 
st=exp {—0.22(A+,,L)} j 


in expressions the same as (8)—(10) but in which the numbers 24, 0.4, 0.2, 
144,°3.35, 2.9, 29, 8:35, 7.9 are replaced by 31, 0:16; 0.251287 3:72.3.1-326- 
8.7, 8.1 respectively, and capital letters are replaced by lower case letters. 

Using the above sets of values, calculations were made using 43 day 
periods. Since T;, (in both cases taken to be 22.5) is so long and since a 
rather large period (43 days) was used in the calculations, instead of 
using the values of the survival functions S7+(t — 5) in (10), the average 
of the values for ¢ — 4, ¢ — 5, and ¢ — 6 was used. Similar averages were 
used for Sz? in (9). In the calculations for T. castaneum the values become 
erratic if this is not done. One other modification was introduced. In order 
not to have to take into account the death rate as it depends upon age, 
but to take into consideration the fact that beginning with the second 
cycle the adults are almost all young, it was assumed that it takes a period 
of 1/D days to reach a homogeneous distribution of ages. Hence the death 
rate was taken to be linear with time beginning at zero when the first new 
adults appeared and being equal to the average value D, after 1/D days. 
Furthermore, adults were assumed to become mature two periods (9 days) 
after emerging, immature adults being assumed to be inactive. The results 
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of the calculations are shown in Figures 1 and 2. It will be recalled that 
the number of pupae and larvae were taken to be represented by the 
number of larvae alone. Hence in the figures the values for the sum are not 
differentiated from those of the larvae alone. In these calculations the 
fractions of the pupae to larvae in the steady state are respectively (Tp = 
6.2, tp = 7.0 days) 0.08/17.6 and 0.45/6.2. 

From an inspection of the figures it can be seen that a fairly satisfactory 
representation is achieved. Note that in the calculated values for T. casta- 
neum the number of larvae falls below the number of adults from about 
90 to 200 days, then exceeds the adult number thereafter. This pattern is 
not so clear in the data, but in the corresponding curves (Park, 1948, ICc, 


OE) 
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Ficure 1. Ordinate: number of individuals per gram of flour, the scale being linear with 
log (1 + NV), where N represents the number of eggs Z, the number of larvae L, or the number 
of adults A. Abscissa: time in days. Curves calculated from equations (8)—(10) together with 
(7). Data on T. castanewm from Park (loc. cit.) ; x—number of adults, o—larvae and pupae. 
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Ficure 2. For description see Figure 1. Curves calculated from equations (8)—(10) together 
with (11). Data on 7. confuswm from Park (loc. cit.). 
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IICc, and IIICc) this is more pronounced. During this period the infection 
is probably just getting started, according to Park, so that this pattern 
should perhaps be the same for all. In the case of T. confuswm it is clear 
that the ratio 1/10 should be replaced by a ratio closer to 1/5 in order that 
the minimum should not be quite so deep. Note, however, in the corre- 
sponding curves for infected populations the values are all lower near the 
minimum. The infection is presumed to have a relatively small effect on 
T. confusum, so the values in these cases should be perhaps about the 
same as in the sterile situation. 

In both cases the values of A,,, . . . /,, were taken from Table 10 of Park 
(loc. cit.). These values are the averages for all times. Because of the tran- 
sient effects it would have been better to use averages in which the first 
three hundred days were omitted. 

It will be noted that the data in Figure 5 of Park’s monograph, of which 
only part is shown in Figure 1, show a tendency for the larval and adult 
curves to oscillate, the two being out of phase. No such effect appears in 
the calculations, except for the crossings at 90 and 200 days. Several pos- 
sibilities for this effect suggest themselves. First, there is the fact that the 
system is not too well damped so that random changes might maintain 
an apparent oscillation. Then there are the effects of age on egg-laying rate 
and death rate which could play a role if the age distribution fails to be- 
come uniform. This would be especially likely if all the adults tend to live 
to be about the same age. The delayed effects of egg-eating on egg-laying 
are only likely to be important if the delay is greater than about a month. 

Since the number of eggs can be calculated readily from expression (8) 
the values are included in Figures 1 and 2 although no experimental values 
were given. In this connection it may be pointed out that from (8) in the 
steady state, but with Rz not replaced by its numerical value, Rz~W/S2? is 
readily estimated and hence using Sz from (5), Rz can be calculated. On 
the other hand, one can also use the egg numbers to adjust Rp and Sz so 
that the prediction of the outcome of competition might be improved by 
using this information in addition to direct estimates of Rp. 

Recently Park (1954) has studied interspecies competition under vari- 
ous conditions. When the values of the various parameters are available 
for each of the experimental conditions, it will be possible to determine the 
extent to which the present model can account for the results of compe- 
tition. 

If three species are mixed in equal proportions then AA is positive if 


Sz SF Se Se eee =) = 1g 
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the o’s being the survival functions for the third species. The condition 
for Aa > 0 is the same except that capitals are changed to lower case, the 
lower case letters are changed to Greek letters and Greek letters are 
changed to capitals, except for the subscripts. 

Competition between two similar species with differential interaction. If 
any stage of one species is more resistant to cannibalism than the corre- 
sponding stage of the competing species, then the survival functions have 
different values for the two species. In general it is necessary to be able to 
measure a coefficient C;; representing the interaction of stage i of species 
2 (lower case) on stage 7 of species 1 (capitals) and a coefficient c;; repre- 
senting the interaction of stage 7 of species 1 on stage 7 of species 2, these 
in addition to the coefficients C;; and c;;. Some estimate of the extent of the 
differences in these coefficients could be obtained from the results of in- 
cubating various combinations of stages from the two species. If there is 
no differential interaction then C,;; = Ci; = ci; = ¢,;. If there is no differ- 
ential resistance to cannibalism then C,; = c;; and Ci; = c,;. In either of 
these special cases the results of the previous section may be used, other- 
wise not. If there is no differential cannibalism, then C,;; = C;; and ¢;; = 
c;;. But these equalities do not help to simplify matters to any extent. To 
predict the outcome of an interaction it is necessary to know each coeffi- 
cient, then substitute the values into the expressions for AA and Aa cor- 
responding to a given mixture after having found the stable state values. 
If, however, it is found that C;;/c;; is the same for each stage i, so that we 
may set these ratios equal to some values C;/c;, we may obtain expressions 
similar to (2) and (3). This case may be referred to as proportional can- 
nibalism. For example, this condition requires that if the adults of a spe- 
cies are twice as cannibalistic toward eggs of their own kind as are the 
larvae, then these adults are twice as cannibalistic as are the larvae to- 
ward eggs of a second species. In this case we may use expressions (2) and 
(3) if we replace Tz by ¢yTgz/cz and Ty; by czT,/cz in the former and & 
by ¢/Tz — ctz, Tz by c’Tz + cfg, etc., in the latter. 

If two species of Tribolium are allowed to interact under various environ- 
mental conditions, expression (3) with the appropriate values of the 
parameters for the given conditions should predict the outcome in each 
case. Furthermore, the model may be useful in the study of other less 
closely related organisms. 

In the above model there has been no consideration of random effects. 
Actually, the variability in rg as well as in the various cannibalistic proc- 
esses would be especially important when the numbers involved are small 
and when the system is weakly damped. Information on these variabilities 
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would, in principle, enable one to estimate probabilities of success in com- 
petition for given values of the index 3. This problem will be discussed in 
a subsequent paper. 


The author is indebted to Drs. Thomas Park and Dennis Strawbridge 
for reading and discussing the manuscript. 

This work was aided by a grant from the Dr. Wallace C. and Clara A. 
Abbott Memorial Fund of the University of Chicago. 
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Many solutions are available to the differential equations for systems consisting of a space 
region with a boundary at which the concentration is fixed, diffusion occurring across this 
boundary. A method is described for readily transforming these solutions into results for similar 
systems in which the diffusing substance is removed by a first-order reaction and also removed 
or produced at a rate which is expressible as a polynomial in the time variable. Subsidiary 
transformations and steady-state conditions are also discussed. An indication is given of bio- 
logical applications of the results made available by this method. 


Processes in which both diffusion and reaction are occurring are funda- 
mental in biological systems. They arise in solute passage between blood 
plasma and erythrocytes, passage through capillary walls into surround- 
ing tissue and into all cells from extracellular fluid. Im vitro and in vivo 
experiments may involve diffusion of substrates into tissues and through 
tissues to sites of enzymic action. Chemical laboratory and industrial 
work may also be affected by these processes, and as the equations of 
diffusion with reaction apply to the conduction of heat with heat produc- 
tion the equations have wide applications in this context. Classical dif- 
fusion only is considered here; the assumption of ideality being most satis- 
factory for small concentrations of diffusing substance in a medium where 
conditions are simple. 

I. General description of the systems. Consider a region possessing just 
one boundary at which net transport of diffusing substances may occur 
and “outside” that boundary the concentration has a fixed value C for 
¢ > 0. Such a situation would be exemplified by the interior of a sphere 
with constant external concentration, or the region exterior to a sphere 
with constant concentration inside, or the annular space between two 
coaxial cylinders with constant concentration within the inner cylinder. 
If there is no resistance to the flux across the boundary there will be a 
surface concentration c* at the “inner” surface of the region, where c* = 
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aC, a being the partition coefficient of the diffusing substance between 
the media inside and external to the region. Thus c* may be taken as the 
effective external concentration and if c’(x, y, 2, 4) = c’(x, ¢) represents 
the concentration in the system then c’ = c* forms the boundary condi- 
tion. If a resistance to the flux, e.g., a membrane, exists then, at the 
boundary, Ddc’/dn = P(c* — c’) where 7 is the normal to the boundary 
directed into the medium of constant concentration, c’ the concentration 
at the surface of the region, D the diffusion coefficient in the region and P 
the permeability of the interface. In accordance with the usage of P. V. 
Danckwerts (1951) and of H. S. Carslaw and J. C. Jaeger (1947b), P/D 
is written as # and the boundary condition becomes 


=h(c*—c). (19 


Another boundary, real or mathematical, to the region will, of course, 
exist but must be such that there is no net flow across it, i.e., at this 
boundary either dc’/dn = 0 or c’ is finite for any finite time. The former 
condition will apply to a definite physical boundary whilst the latter is 
most convenient for such boundaries as the central point of a spherical 
region, the axis of a cylindrical region, or at infinity. 

It is assumed that a uniform concentration co exists initially within 
the region. The literature and textbooks (e.g., Carslaw and Jaeger, 1947a) 
contain solutions, giving c’ as a function of space coordinates and time 
for regions of many shapes, particularly in the special cases, included in 
all the present discussions, where cy = 0. Under the latter circumstances 
the solutions for c’ are given the symbol y as discussed in Section II be- 
low. Results for c’ for any initial concentration cy can always be obtained 
from the corresponding solutions y, simply by adding co to the function 
obtained by replacing c* by (c* — co) in the expression for y. 

IT. Transformations. In this section several expressions are commented 
upon without giving proofs. The proof of the most important of these ex- 
pressions is given in the following section and the other results mentioned 
may be verified in a similar way. First consider the transformation 


c= c' exp(— kl) +h f 0x, t) exp (— kt) di 


1 x t t 
ee myn\ f LG, t) exp (— kt) (di)*+1 Q) 


+ exp (— ki) a et man kata 
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where & is a constant and the m,’s are constant coefficients. The symbol 
refers to the known solution c’ to one of the above diffusion problems in 
which zero is substituted for co, if this is not already zero. The integrands 
are strictly to be considered as functions of a variable of integration ¢’ 
rather than of ¢ itself. Application of this transformation gives the solu- 
tion c to the similar problem in which not only the diffusion process oc- 
curs but also, in the region, the solute is removed by a first-order reaction 
(rate kc) and is produced or removed at a rate which is expressible as a 
polynomial function of time. A very wide variety of time variable rela- 
tions may be, at least approximately, represented in this way. 
Thus from c’ the solution to the equation 
dc’ 
ot 


= DY’ c’, (3) 


with the initial condition c’ = cy at t = 0 (c’ = 0 being included) and the 
following conditions: either c’ = c* or dc’/dn = h(c* — c’) at one bound- 
ary, and either dc’/dn = 0 or c’ finite for finite time at the other bound- 
ary; the transformation gives c, the solution to 


N 
Gr = DUC hot Se mal (4) 
for the same initial and boundary conditions in c. 

The function y, obtained from c’ simply by giving the value zero to co 
in the mathematical expression for c’ in any particular case, may have its 
precise relation to c’ illustrated as follows. For example, if c’(x, ¢) is a 
solution of equation (3) with 


etr, ON= co = const. 


So ]ene Let xan, 


ey =e 
Onix=x2 


where x, and x; are points of the first and second boundary respectively, 
then y/(x, ¢) will be a solution of 


= DVY with w(x, 0) =0, 
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Epes 


If & is zero, the transformation 


, 1 a : j n+1 SS tint 5 
cmc mae marl fs f ve) (a0 Sheen (5) 


applies and may be treated as for the most general case. Where it appears 
in applications, & is usually positive in sign, but equation (2) does not 
carry this restriction. Many reactions which are not first order may be 
approximately represented in this way under suitable conditions (Danck- 
werts, 1950) and if a reversible first-order reaction is involved the results 
will apply to the early stages of the process (e.g., Roughton, 1952a). 
Cases of diffusion problems with reversible first-order removal reactions 
have been investigated by J. Crank (1952). The polynomial may be 
merely a positive or negative constant or may not be present at all. In 
the latter case, all the coefficients m, are zero and the right-hand side of 
equation (2) simplifies to just the first two terms (cf. Danckwerts, 1951). 
If the polynomial is ever negative, care must be taken in that, whilst 
negative values for c are permissible mathematically and will have sig- 
nificance for conduction of heat problems where an arbitrary zero of 
temperature has been taken, they have no significance for concentrations 
and, in practice, the removal operation will not occur if at any place or 
time c becomes zero. 

For those cases where c* has the value zero, equation (2) may some- 
times usefully be replaced by 


and 


ee t t 
was / ames ee , ee. n 
c=c exp(—hki) + ape mant fe a fh c exp (— ki) (di) ** 


which is equivalent to (2) only in these cases. The numbers & and m, may 
adopt any value including zero and the transformation (2) reduces to 
c = c’ exp (—At) if all the numbers m, are zero and only if c* is zero. 
If, under the conditions considered in this paper, c’ is the solution to 
oe, 
Ot 


= DV? c! + pon : 


where » is a positive improper fraction and m the greatest whole number 
contained in v, then equations (2) and (5) respectively provide the solu- 
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tions to 


and to the corresponding equation where k = 0. 
III. Proof that equation (2) is a valid transformation. The operation 
DY’ applied to equation (2) gives 


t 
DY?c =exp (— kt) DV? +k f exp (— kt) DV*ydt 
0 


1 N t t 
— ! po n 
Fee mont ff sea f exp ( kt) DV2y (dt) =*1 


and on substituting for Dy’c’ from equation (3), and similarly for Dy’ 
this becomes 


ae UR Ae 
DY? c =< exp ( kt) +h f exp (— ki) di 


1 ¥ t oy 
=e = n+1 
ee > myn J cee fl i exp ( kt) (dt) 5 


By differentiating equation (2) with respect to ¢, adding kc obtained from 
equation (2), and suitably rearranging the terms we get eventually 


of+k cae exp (—kt) +2[¥ (x, 4) exp (—R?) +h fy (x, t) exp (— kt) at| 


> mart fo [9 Dex (— kD a+. {YH 


N n m—r—1Bb—r— Oy Oe 
Xexp (—k?) heaps mnt! > (1) [7+ = 7) iE 


where the quantity 


fo ti 


(n—r—1)! 


must be given the value zero for r = 1 since this term is obtained as the 
derivative of a constant. The law for integration by parts can then be 
applied to the first two square bracketed terms. Also the inner sum of the 
final double sum is independent of r and by writing out the first few terms 
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is seen to be #"/m!. Thus dc/dt + ke becomes 


Oo exp (= bi) +f oY exp (— Bi) dl 


-2> m nf. Lf exp (- 8) (att SS my 
9 n=0 0 dt n=0 


and from equation (6) this is seen to be 


DV? c+ 3) Mt”. 


n=0 


So if c’ is a solution to equation (3), then c, defined by transformation 
(2), is a solution to equation (4). 


Initial condition. If c’ = co at t = 0 then from equation (2) 


N M,N! N = ( = 1) tin pai 
BD re Rye t lim ne mint! De (n = r) ! 


Boundary conditions. For the consideration of these the following result 
is required 


fo. fep(-b) dns op + ~~ bee) Deeley age 


(a — ta" al 


This is most simply proved by induction. Direct substitution shows that it 
holds for n = 0 and m = 1. Taking m = s and integrating both sides of 
equation (7) with respect to ¢ gives 


ff el (dt) s+2 


_exp(— Ai) | (—1) 1) SK (=) rong 
net et ria 


r=0 


exp Ce at) st (—1) "isti-rp-r-1 
paleak ete a (sti—vr)h” 


Thus if equation (7) holds for » = s, it is also true form = 5 + 1, which 
establishes the required result. 
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Now if c’ = c* at the boundary, then 


t 
c= c* exp (— kt) +ho* f exp (— kt) di 
0 


=> ment fi. ; ff exp(—#i) (dt) +1 
n=0 


' 2s rjn—r pb—r— 
+ exp (— ki) SS ( DDS oo : 


By applying result (7) and evaluating the first integral we find that c is 
also equal to c*. | 

Taking the directional derivative of c with respect to the normal 7 to 
the boundary surface gives 


fe) 
$e = oo exp (= Bi) thf s se exp ( kt) di 
-4> y marl fo Ape se exp (— Bi) (di)™*1, (8) 


If equation (1) represents the boundary condition on c’, then on substitut- 
ing for dc’/dn and dW/d7n in equation (8) and rearranging one obtains 


oo = ho — hl o! exp(— kd) +k f v(x, t) exp (— ki) dé 
0 


on 
pea t t 
aie myn f veel, i) exp (— kf) (dt)"#1 


+S mnt fi.  f exp (— 2d) (deat], 
n=0 


On applying result (7) to the last summation it follows that 


i =h(c*—c). 

For conditions at the other boundary, it is clear from equation (8) 
that where 0c’/dn (including dW/0n) is zero, then 0c/dy is also zero. Final- 
ly, it follows from the original transformation that, if c’ is finite at a 
boundary, then c is also finite at this boundary. 

IV. Non-commutative nature of the transformations. Transformation (2) 
may be effected in stages in a definite way. For the simplest and most 
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useful case, consider solutions, under our specified conditions, to the 
equations in the following scheme, starting with equation (3): 


ae! 
ot 
a — 
OC 


Pap ee mo 2 DY kG (9) 


Soles 


—ke + Mo. (1 0) 
By appropriate application of equation (5), transformation 7, may be 
effected, i.e., c: may be obtained from c’. Now it can easily be shown that 
c is obtained by direct application of equation (2) with all m, values 
zero, i.e., transformation J», to all the function ¢,. The application of T> 
followed by 7; does not transform c’ into c; this process may, however, be 
effected by following 72 with a transformation J” defined by 


=s 


t 
c= at f c (¢* = 0) dt, 


where ¢2(c* = 0) represents the solution c, in which c* is given the value 
zero. These results may be expressed symbolically in the form, 


T, eli — tia XL Ls <n 5 


V. Rate of uptake of diffusing substance. For many problems, e.g., dif- 
fusion into cells, the rate of uptake and the quantity of substance ab- 
sorbed are of interest. If no chemical reaction is involved the rate of up- 


take will be 
RAG = f [05 as=5 ff f cev, (11) 


Q’ being the quantity which has been absorbed in a given time. When re- 


actions are in progress, 
r= = [fo yas, (12) 


where either R or Q can, of course, be negative. 
If R’ is known, then R may be found in any of the problems considered 
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in this paper. For example, from equations (2), (11), and (12) 


t 
R=R' exp(— kl) +k f Wexp(—kd) di 
0 


BLY ae man fr. ale W exp (— ki) (di)7*!. 

The function W is the rate of uptake in the nonreacting system when 
the initial concentration is zero and it has a similar relation to R’ as y 
has to c’. 

VI. Steady-state theory. In many diffusion systems, as ¢ increases, the 
concentration tends toward a constant value at each point. The steady- 
state solutions may be obtained by letting ¢ tend to infinity in the solu- 
tions for the transient state. This procedure, however, does not usually 
give the results in the most convenient form. Now, under the conditions 
considered in this paper, a steady state will exist if equations (9) or (10) 
apply, with & > 0 in each case. An ingenious method for obtaining solu- 
tions in the former case has been discussed by Danckwerts (1951). His 
results obtained for the initial condition co = 0 apply without modifica- 
tion if an initial concentration is present. 

The transformation T, obtained from equation (2), giving the solution 
to equation (10) is 


¢ = c’ exp(—&i) + (4-2) fv, t) exp (— ki) di 


+=*{1—exp(—#i) }. 


If ¢ tends to infinity then 
oe m 
to = (2-72) f v(x, Dexp(— kd) dit. 


The integral, assuming it to exist, is the Laplace transform (x, &) of 
the function (x, #), considered as a function of k instead of p (Carslaw 
and Jaeger, 1947e). Thus the steady-state concentration is given by 


coo = + (bE) Ox, B) (13) 
and, similarly, the rate of uptake is 
Ro = (#2) ¥(h) (14) 


Thus, in the cases discussed, if the required Laplace transforms happen 
to be to hand, the steady-state solutions may be written down immedi- 
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ately. If the transforms have to be calculated, this necessitates solving 
a differential equation and no advantage is gained over solving the original 
differential equation with 

Oc 


ay 


VII. Application of method with examples. Solutions to the equations of 
classical diffusion processes are usually obtained by the method of 
Laplace transforms (Carslaw and Jaeger, 1947e). The application of this 
method is elementary and rapid when the transforms involved occur in 
tables, the required results being usually given in terms of tabulated 
functions. Otherwise the inversion theorem has to be applied and this 
requires some skill at contour integration and may involve much compu- 
tation before results, usually expressed in terms of improper integrals or 
infinite series, are obtained. Under the latter circumstances identical re- 
sults are given by the easily applied transformations described in this 
paper whilst in the former the results are not usually expressed directly 
in terms of tabulated functions. Transformations of the type considered 
here can also give results in cases where difficulties arise in the Laplace 
transform method. This is to be discussed in a later paper. 

With the present problems, the function (x, #) frequently takes the 
form c* — c*d, where ¢ is a function of the space coordinates and time. 
This can result in considerable simplification as equation (2) then be- 
comes 


t 
c= c! exp(— ki) + o*}1—exp(— kd) —k f sexp(—kt att 
0 


N t t 
+ Do inn f fb exp (— ki) Oe. 


where the numbers & and m, can have any values, including zero. Usually 
¢ is either an infinite series of terms or an improper integral of a function, 
in which the time variable occurs only in a factor of type exp (—A/). 
Thus integrations with respect to time of ¢ exp (—At) are readily per- 
formed. 

Two examples are now given of the application of these transformations. 

(a) Diffusion in a direction normal to the faces of an infinite plane slab 
(—1<x <]) with zero initial concentration and a constant concentration c* 
atx = +1. If no reaction occurs then 


Aa Se cos x (sy sexp (— ut) , 


(15) 


c=(x,) = ch 
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where u = (2s + 1)?1D/4/? (Carslaw and Jaeger, 1947c). 
For the similar system with production of diffusing substance in the 
slab at rate mo + mut the required transformation, from equation (15), is 


c= elt me f adit fo fo (dn? 


Application of this gives 


— 
where 


m m 
v= (wer + m—™) exp ( — u1) =i or mt. 


Now R’ (the rate of uptake through unit area of one face) is here Ddc’/dx 
at « = / and this becomes 


4¢* Di2— (—1) 21/2 : (Sp ie 
r ~~ Peet sin |! D {+ exp ( ue 


8=0 


Application of the transformation 


t t t 
R=R-2f Ral f Ran? 
0 


gives 
ei De— ( == 1) \1(* yt. 
eT: ores ace NAGS 


Tf desired, Q = f Rdi may then, immediately, be obtained. 

(b) Diffusion into a sphere (0 <r < a) with an initial concentration co 
and a resistance to the flux at the surface of the sphere. If no reaction occurs, 
c’ is given by the solution of the corresponding problem where initially 
c’ = 0 (Danckwerts, 1951) plus the solution when c’ = ¢p initially and 
c* = 0 (Carslaw and Jaeger, 1947d), ie., 


co 
ot = ot — 2 (cs — 6) »> 2, exp (— Dai?) «sin aa, -sin Ta 
s=1 a; 


where 
ie a’a2+ (ah—1)? 
aa oa eee: 
B=ae2+h(ah—1), 
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and a, are the successive roots of the equation aa cot aa + ah —1=0 


(Carslaw and Jaeger, 1947f). 
For the solution to the corresponding problem in which first-order 


removal together with zero-order removal (—ke — m) occurs in the 
spherical region, the transformation 


t 
c= cl exp(— kt) + c*{1—exp (— kD} — (ko*+m) f ¢ exp(— kf) dl 
0 
from equation (15), is required. This gives 


2h ss 9 


y Kamd 0,2 6 
s8=l 8 


Gas 


[ (c* Da? =n —¢; bexp (01) ke a 


— 
Xsin @a,-sin ras , 


where 6 = Da? + k. The rate of uptake by the sphere when no reaction 
occurs is 


/ 
Roa (<5) 
Olena 


= 81a?Dh? (c*— co) >) a exp(— Da) , 
3s=1 


RlrR 


and, applying the transformation, 
t 
R=R’ exp(— kt) +(i +) if W exp (— ki) dt 
c 0 


gives 


| 


R= 8nra?Dh? >) —[(c* Da2— m— c,5) exp(— 6f) +ke*+m] . 
s=1 


DR 


6 


Also 0 becomes 
8ra°Di > 1. (het + m) dt-+ (c* Dat — 8) {1 
a BP a2 — m— 6b) (1 —exp(— ede 


Transformations (13) and (14) will apply to the steady state which 
arises here. Thus from the Laplace transforms we get 


kad. ote (kc*+ m) ha? sinh (rp) 
rk{aucoshaun+ (ah —1) sinh ap} 
and 


ee 4ra°hD(kc*+ m) (au cosh au —sinh ap) 
k{ au cosh au+ (ah —1) sinh ap} i 
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where up = (k/D)"?. Steady-state diffusion equations for spherical! celis 
have been solved by N. Rashevsky and approximately by H. D. Landahl 
(cf. Rashevsky, 1948). 

This system forms a model for the diffusion of oxygen into the spherical 
form of the red corpuscle with first-order reaction with hemoglobin, as- 
sumed to be uniformly distributed in the cell, and allowing for some 
direct metabolic utilization of the oxygen within the cell. The plasma 
membrane is assumed to act as a simple surface resistance and to possess 
no thickness. Carbon monoxide uptake is also covered if m is taken to 
be zero. The equations given here might more accurately describe this 
system than those of F. J. W. Roughton (1952b) as the permeability of 
the membrane has not been neglected. 


The author wishes to thank Professor Sir Charles Dodds, F.R.S. for 
the provision of certain facilities. 
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